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EXECUTIVE  SUMMARY 


Sailor’s  thermo-physiological  responses  to  the  unrelenting  heat  accumulation  possibilities  of  a 
disabled  submarine  were  simulated  with  USARIEM  thermoregulatory  models.  The  objective  of 
this  analysis  was  to  estimate  the  probable  time  and  interior  environmental  conditions  (air 
temperature  and  humidity)  when  heat  related  illnesses  begin  to  occur,  and  the  time  available  for 
successful  rescue  of  the  individuals  after  they  put  on  escape  suits.  The  modeling  indicates  a 
sedentary  lightly  clad  sailor,  waiting  for  rescue,  would  begin  to  experience  thermo-regulatory 
failure,  i.e.,  a  rapid  rise  in  body  temperature,  at  an  air  temperature  (Ta)  about  92°F  (33°C) 
together  with  high  humidity  (90-100%  RH).  From  that  point  on,  the  lightly  clad  sedentary 
individuals  can  survive  for  2,  3  or  6  more  days,  as  air  temperature  continues  to  rise  at  rates  of  3, 
2  or  1  °F  /day  (1.7, 1.1  or  0.6°C/day)  respectively,  before  the  individual’s  body  core  temperature 
(Tc)  reaches  the  unsafe  103°F  (39.5°C)  level.  Donning  an  impermeable  escape  suit  eliminates 
the  evaporative  cooling  from  sweating,  causing  Tc  to  rise  even  more  rapidly.  At  92°F  (33°C) 
conditions  and  wearing  the  escape  suit,  Tc  will  reach  103°F  (39.5°C)  in  about  2  hours.  If  the  suit 
is  put  on  earlier  at  a  cooler  interior  temperature  of  86°F  (30°C)  the  time  available  before 
reaching  103°F  (39.5°C)  increases  to  1  -  3  days  depending  on  the  rate  of  rise  in  ambient  air 
temperature.  Further,  knowing  cabin  air  temperature  and  humidity  enables  estimates  to  be 
made:l)  for  the  beginning  of  Tc  thermoregulatory  failure,  and  2)  of  Tc  levels  of  101. 4°F 
[38.5°C],  103°F  [39.5°C],  and  105°F  [40.5°C])  (Figure  15). 
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INTRODUCTION 


U.S  Navy  submarines  rarely  become  disabled.  The  first  disablement  and  loss  occurred  in 
1914  when  leaking  battery  acid  in  USS  F-4  corroded  hull  rivets  that  caused  progressive  flooding, 
chlorine  off  gassing,  loss  of  depth  control  and  eventual  catastrophic  hull  failure.  Then,  after 
another  accident  in  1927  when  the  U.S.  Submarine  S-4  became  disabled  and  was  lost  with  all 
hands,  rescue  procedure  planning,  equipment  and  training  have  continuously  been  undergoing 
improvement.  As  a  result,  12  years  later,  in  1939  off  the  shore  of  Portsmouth,  N.H.  Sailors  were 
successfully  rescued  from  the  USS  Squalus,  which  quickly  sank  after  a  rear  compartment 
flooded  (Stewart,  2008).  The  disabled  Squalus  was  located  on  the  sea  floor  at  a  depth  of  240  ft  in 
29°F  (-1.7°C)  water,  a  rescue  ship  with  a  diving  chamber  came  to  the  site  and  the  33  crew  in  the 
non-flooded  compartments  were  transferred  to  the  surface  within  40  hours  via  4  trips  of  the 
diving  chamber.  A  major  thermal  hazard  in  such  a  disabling  has  been  the  threat  of  hypothermia 
from  the  environment  caused  by  cold  ocean  water  while  waiting  the  sub’s  rescue  (Castellani  et 
al„  2005). 

In  recent  years,  the  operations  in  warm  sea  water  locations  and  evolving  submarine  designs 
have  increased  the  possibility  of  rising  interior  air  temperatures  in  a  disabled  submarine  with  the 
threat  of  sailors  developing  hyperthermia  and  related  heat  illnesses.  Simulated  disabling  tests 
performed  with  the  USS  Dallas  near  Groton,  CT  in  2004  and  USS  Salt  Lake  City  off  the  San 
Diego  coast  in  2006  found  that  the  interior  air  temperature  and  humidity  slowly  but  consistently 
increased  (Horn  and  Daniel,  2007).  While  disabled,  the  USS  Salt  Lake  City’s  interior 
temperature  rose  from  72  to  85°F  (22  to  29°C)  and  the  humidity  rose  from  50  to  85%  RH  over  4 
days. 

In  the  present  analysis,  sailor’s  thermo-physiological  responses  were  simulated  with 
USARIEM  thermo-physiological  models  to  predict  a  sailor’s  health  risks  from  slowly  rising 
temperature  and  humidity  that  may  occur  in  a  disabled  submarine.  The  simulated  predictions  of 
the  sailors’  physiological  responses  could  then  help  guide  rescue  and  planning  operations. 


METHODS 


Modeling 

The  principle  model  used  in  the  present  analysis  was  SCENARIO  (Kraning  and  Gonzalez, 
1997)  because  it  has  the  most  appropriate  physiological  mechanism  features  for  this  application. 
Also,  another  related  thermoregulatory  model,  the  Initial  Capability  Decision  Aid  (ICDA) 
(Yokota  and  Berglund,  2006),  was  customized  to  estimate  the  accumulation  of  water  vapor  in  the 
air  of  a  disabled  submarine  during  rising  temperature  situations  and  correspond  to  that  measured 
in  the  tests  of  the  disabled  USS  Salt  Lake  City. 

SCENARIOS 

SCENARIO_J  is  the  Java  language  version  of  the  SCENARIO  thermo-physiological  model 
developed  at  USARIEM  by  Kraning  and  Gonzalez  (1997)  (Fig  1).  It  is  a  rational  human 
thermoregulation  model  developed  from  first  principles  of  physiology,  heat  transfer  and 
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thermodynamics.  The  physiology  and  heat  transfer  principals  and  properties  are  based  on 
extensive  laboratory  experiments  from  the  literature  and  USARIEM.  The  main  mechanisms  for 
active  physiological  body  temperature  control  are  regulating  sweating,  shivering  and  blood  flow 
to  the  skin.  The  model  accounts  for  the  anthropometries  of  height,  weight  and  %body  fat.  It  also 
accounts  for  the  effects  of  dehydration  and  heat  acclimation. 

Figure  1 .  Schematics  of  SCENARIO_J  model  with  compartment  arrangements  and  weights  for  a 
standard  man  (body  weight,  70kg;  body  surface  area,  1.8  m2). 


The  model  represents  the  human  as  six  lumped  parameter  tissue  compartments  (Figure  1), 
five  of  which  are  concentric  cylindrical  layers  representing  core,  muscle,  fat,  vascular  skin  and 
outer  skin.  The  sixth  is  the  central  blood  compartment  which  interacts  and  exchanges  heat 
directly  with  all  of  the  compartments  except  the  outer  skin.  At  rest,  the  core  is  the  principle 
source  of  heat.  For  activity  levels  above  rest,  the  increased  energy  is  generated  in  the  muscle 
compartment.  With  exercise,  heat  from  working  muscles  predominates  and  can  be  more  than  10 
times  that  of  the  core.  Heat  flows  by  conduction  between  adjoining  compartments  and  by 
convection  through  blood  flow  from  the  central  blood  compartment.  These  heat  flow  pathways 
help  to  distribute  the  heat  to  the  other  compartments.  The  environment  exchanges  heat  and 
moisture  primarily  with  the  skin  (through  clothing  and  exposed  skin)  and  also  directly  with  the 
core  through  respiration. 

SUB_ICDA 

The  Sub_ICDA  model  is  a  version  of  the  USARIEM  ICDA  model  (Yokota  and  Berglund, 
2006)  customized  for  this  submarine  application.  ICDA  is  an  adaptation  of  Gagge  (1971,  1986) 
model  and  SCENARIO_J.  The  ICDA  model  is  similar  to  SCENARIO_J  but  simpler  with  the 
human  represented  as  two  physiological  compartments  (core  and  skin). 

All  of  the  model’s  metabolic  activity  (M)  takes  place  in  the  core  compartment  (Figure  2). 
The  customization  enables  inputting  metabolism  directly  using  the  met  unit  to  characterize  and 
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quantify  metabolism  (met  =  actual  metabolism/resting  metabolism).  The  model  then  estimates 
the  heart  rate  necessary  for  the  metabolism  and  thermoregulatory  cardiovascular  responses 
related  to  the  environment.  The  sailor’s  evaporative  and  respiratory  water  losses  are  used  to 
estimate  humidity  level  changes  in  the  occupied  submarine  environment.  Further  Sub_ICDA 
details  are  available  in  the  Appendix  A. 

Figure  2.  Schematic  of  Sub_ICDA. 


=  core 

I  l  =  skin  compartment 
I — l  =  clothing  property 


Respiratoryheat  loss 


Evaporation  from  skin 


Skin  blood  flow 


Convection  to  air 


Conduction  through  clothing 


Radiation  to  surroundingsurfaces 


SUBMARINE  ENVIRONMENT 

The  environment  of  occupied  spaces  in  a  disabled  submarine  can  change  in  different  ways 
depending  on  the  situation  and  cause  of  the  disablement.  This  project  considers  the  submarine  to 
be  powerless  but  intact  and  unable  to  control  temperature,  humidity  and  air  motion  after  being 
disabled.  Further  the  project  considers  the  temperature  of  the  interior  surfaces  to  approximately 
equal  the  occupied  space’s  air  temperature.  The  disablement  details  are  not  considered  but  only 
the  human  response  to  slowly  rising  temperatures  that  follow  the  incident.  Also,  the  primary 
source  of  water  vapor  to  increase  air  humidity  is  considered  to  be  that  from  the  occupants.  The 
water  vapor  generated  by  CO2  scrubbers  adds  about  16  g/h  per  resting  sailor  and  was  neglected. 

The  disablement  test  of  the  USS  Salt  Lake  City  measured  the  interior’s  temperature  and 
humidity  for  a  period  of  four  days  following  the  disabling.  The  temperature  and  relative 
humidity  increased  rather  linearly  from  76°F  (24.4°C)  and  60%  RH  to  84°F  (29°C)  and  85%RH 
during  the  4  day  test  interval.  An  earlier  disabling  test  on  the  USS  Dallas  found  a  similar 
temperature  rise  over  a  3  day  interval,  though  the  humidity  rise  was  slower  and  less  linear 
(Figure  3)  (Horn  and  Daniel,  2007). 
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The  SCENARIO_J  model  needs  temperature  and  humidity  as  inputs.  It  is  assumed  that  the 
temperature  in  disabled  subs  will  rise  approximately  linearly  as  the  test  measurements  of  Figure 
3  indicate.  Further,  it  is  likely  that  the  temperature  will  eventually  become  asymptotic  to  some 
upper  limit  when  the  vehicle  reaches  a  quasi-steady  energy  balance  state  depending  on  the 
conditions  of  the  sub,  its  crew  and  the  temperature  of  the  surrounding  seawater.  At  present,  the 
upper  level  temperature  limit  mechanisms  are  undefined  and  unpredictable.  Instead,  the  interior 
temperature  is  assumed  to  rise 

Figure  3.  Measured  temperature  and  humidity  of  occupied  spaces  in  disabled  submarines. 
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Figure  from  Horn  and  Daniel,  (2007) 

steadily  without  an  upper  limit  from  the  pre-disablement  temperature.  The  humidity  may  also 
rise.  Humidity  was  predicted  with  the  Sub_ICDA  model  for  different  rates  of  temperature  rise, 
assuming  the  water  vapor  in  the  air  increases  only  from  respiration,  evaporated  sweat  and  water 
diffusion  through  the  skin.  In  turn,  water  in  the  air  can  be  absorbed  by  fabrics,  wood  and  various 
surface  materials.  To  model  the  humidity  rise  for  various  situations,  an  effective  air  volume 
(Veff)  or  air  mass  per  person  was  determined  by  trial  and  error  to  correspond  to  the  humidity  rise 
observed  in  the  disabled  USS  Salt  Fake  City.  This  effective  air  volume  could  then  be  used  to 
predict  the  humidity  behavior  for  other  rates  of  temperature  rise  and  starting  points.  Some  of  the 
trial  and  error  results  are  displayed  in  Figure  4  for  the  2°F  (1.11  °C)/day  rise  from  76°F  (24.4°C) 
and  60  %  RH  to  find  the  Veff  that  results  in  model  predicted  RH  levels  that  correspond  with 
those  measured  in  the  disabled  USS  Salt  Fake  City. 
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This  effective  volume  of  360  m  /person  was  used  by  Sub_ICDA  to  determine  the  humidity 
behavior  for  other  rates  of  temperature  rise  and  starting  points.  The  effective  volume  is  larger 
than  the  actual  air  volume  per  sailor  because  it  agrees  with  the  disabled  submarine  test 
measurements  and  thereby  also  accounts  for  unknown  adsorption  and  condensation  mechanisms. 
Results  for  rates  of  1,  2  and  3  °F  /day  (0.56,  1.11,  and  1.67  °C/day)  starting  from  76°F  (24.4°C) 
are  plotted  in  Figures  5  a  and  5b.  The  relative  humidity  responses  with  time  (Figure  5  a)  are  very 
similar  for  the  1,  2  and  3°F/day  rates  of  air  temperature  rise;  reaching  100%  in  about  6  days. 
However  the  absolute  humidity  (Figure  5a)  or  water  vapor  concentration  levels  in  the  air 
(lb_w/lb_dryair)  continues  to  rise  after  the  relative  humidity  reaches  100%  because  the 
temperature  continues  to  increase,  allowing  the  air  to  absorb  more  water  vapor.  The  relative 
humidity  and  time  results  of  Figure  5a  are  used  by  SCENARIO_J  to  predict  the  thermoregulatory 
responses  of  sailors  to  these  conditions.  Figure  5b  indicates  how  the  dew  point  (Tdp)  of  water 
vapor  in  the  air  increases  until  it  reaches  air  temperature  and  then  increases  with  Ta. 

Figure  4.  Results  of  trial  and  error  method  to  estimate  effective  air  volume  (Veff )  per  person. 

Simulated  (solid  lines)  and  measured  (dotted  line)  relative  humidity  in  disabled  submarine  USS 
Salt  Lake  City  with  2°F  (1.2°C)/day  temperature  rise  from  76°F(24.4°C). 


Figure  5a.  The  relative  humidity(RH)  and  absolute  humidity  or  water  concentration  in  the  air 
generated  by  sailors(Hughes  and  Horn,  2006)  resting  in  shorts  and  T-shirts  in  a  disabled 
ubmarine  with  air  temperature  increasing  at  different  rates  starting  from  76°F 
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Figure  5b.  Dew  point  temperatures  (Tdp)  resulting  from  sailors  resting  in  shorts  and  T-shirts  in  a 
disabled  submarine  with  air  temperature  (Ta)  increasing  at  different  rates  starting  from  76°F 

(24.4°C)  and  60%  RH. 


PHYSIOLOGICAL  ASSESSMENT 

Core  temperature  (Tc),  skin  temperature  (Tsk),  sweat  rate  (SR),  and  heart  rate  (HR),  are 
predicted  by  the  models  in  response  to  the  environment,  clothing,  activity  and  time  and  indicate 
the  physiological  state  of  the  individual.  In  addition,  the  physiological  strain  index  (PSI)  (Moran 
et  al.,  1998)  is  also  evaluated  to  assess  thermal  and  cardiovascular  strain  and  the  potential  for  heat 
illness  based  on  the  ratio  of  deviations  in  Tc  and  HR  from  resting  thermally  neutral  levels  to  the 
maximum  deviations  acceptable  for  the  individual,  where: 


PSI  =  5  [(HR-HRrest)/ (HRmax-HRrest)]  +5[(Tc-Tc_rest)/(Tcmax-Tc_rest)] 


HRmax  is  the  maximum  heart  rate  for  the  person’s  age  (HRmax  =  220  -  age);  Tcmax  =  39.5°C,  used 
in  these  simulations,  is  the  level  where  the  incidence  of  heat  illness  is  likely  to  occur  in  50%  or 
more  of  males  (Sawka  et  al.,  2001);  Tc_rest  =  36.9°C,  corresponds  to  the  core  temperature  of 
persons  at  rest  in  a  thermally  neutral  environment  that  requires  no  thermoregulatory  activity. 
HRrest  =70  bpm  is  the  rate  at  rest  in  a  thermally  neutral  environment  used  in  these  simulations, 
together  with  an  age  of  23  yr.  Thus  PSI  is  10  when  Tc  and  HR  equal  maximum  levels.  The 
physiological  strain  scale  between  0  and  10  generally  corresponds  to  a  categorical  scale  from: 
None/slight  to  Low,  Moderate,  High,  and  Very  High  as  in  Fig  6f.  PSI  levels  above  8  represent 
very  high  strain  states  and  are  to  be  avoided  except  for  short  intervals. 


HEAT  ACCLIMATION 

The  effect  of  heat  acclimation  on  sailor  responses  to  conditions  (Fig  6a)  in  the  disabled  USS 
Salt  Lake  City  as  estimated  by  SCENARIO_J,  are  displayed  in  Figures  6b,  c,  d,  e  and  f  The 
graphs  indicate  that  full  heat  acclimation  (fac)  is  beneficial  in  reducing  heat  strain  initially. 
However,  in  the  region  where  core  temperature  (Fig  6b)  starts  to  rise  rapidly  as  thermoregulatory 
failure  begins  and  beyond,  there  is  little  difference  in  modeled  responses  between  acclimated  and 
un-acclimated  (uac)  sailors. 
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Figure  6a.  Ambient  conditions  in  submarine  modeled  after  measuremets  from  disabled  USS  Salt 

Lake  City. 
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Figure  6b.  Core  temperature  response  to  conditions  of  Fig  6a  for  resting  heat  un-  acclimated 
(uac)  and  acclimated  (fac)  50th  percentile  sailors. 


Figure  6c.  Skin  temperature  response  to  conditions  of  Fig  6a. 
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Figure  6d.  Heart  rate  response  to  conditions  of  Fig  6a. 


Figure  6e.  Sweat  rate  responses  of  resting  un-acclimated(uac)  and  fully  acclimated  (fac)  50th 
percentile  male  sailors  during  conditions  of  Fig  6a. 


Figure  6f.  Heat  acclimation  effects  on  the  physiological  strain  index  for  resting  un-  acclimated 
(uac)  and  fully  acclimated  (fac)  sailors  during  2°F/day  rise  conditions  (Fig  6a). 
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INDIVIDUAL  VARIABILITY  -  Body  Size  Effects 

The  U.S.  Navy  consists  of  individuals  with  diverse  body  sizes  and  composition  (Table  1) 
(Hughes  and  Horn,  2006)  (Yokota,  2011)  and,  consequently,  the  physiological  responses  of 
individuals  within  the  Navy  populations  may  differ  from  average  responses.  Assessing  the 
thermal  responses  of  different  individuals  is  important  for  rescue  purposes  and  to  create  safety 
guidelines  for  military  populations.  Submariners  live  for  long  periods  in  constant  comfortable 
conditions  and  are  unlikely  to  be  heat  acclimated  or  to  retain  it.  Further  as  figures  6b-f  indicate 
acclimatization  has  little  effect  once  thermoregulation  begins  to  fail.  Thus,  this  study  examined 
the  effects  of  body  size  and  composition  on  un-acclimated  (uac)  sailors. 

Table  1.  US  Navy  Percentile  Anthropometric  Summary 


percentile 

5th 

50th 

95th 

Height  -  inch  (cm) 

64.7  (164.3) 

69.4  (176.3) 

74.6  (190) 

Weight  -  lb  (kg) 

133.4  (60.5) 

178.6  (81) 

227.3  (103.1) 

Body  fat  -  % 

9.81 

16.9 

24 

The  effects  of  these  anthropometries  on  thermo-physiological  responses  were  estimated  by 
SCENARIO_J  for  the  5th,  50th  and  95th  percentile  sailor  experiencing  the  conditions  (Fig  6a)  of 
the  disabled  USS  Salt  Lake  City.  The  results  are  displayed  in  Figures  7a-e. 


The  core  temperature  (Fig  7a)  and  PSI  (Fig  7e)  responses  indicate  that  the  smaller  (5th  percentile) 
sailors  can  maintain  a  reasonable  steady  body  temperature  longer  in  the  rising  heat  than  the 
bigger  (95th  percentile)  sailors.  The  difference  is  about  a  day  for  a  compartment  temperature 
rising  at  2°F/day  which  may  be  useful  in  designating  the  order  of  escape  to  reduce  possible  heat 
illness  and  its  effects  on  egress.  The  results  also  illustrates  that  individual  water  loss  (sweat  rate, 
Fig  7d),  and  the  need  to  replace  it,  progresses  slowly  until  about  the  8th  day  and  the  beginning  of 
a  rapid  increase  in  body  temperature  when  Ta  is  about  33°C  (91°F).  This  sweat  rate  versus  time 
or  air  temperature  response  patterns  may  suggest  the  rate  of  water  distribution  to  the  crew  if  the 
water  supply  is  limited. 

Figure  7a.  Core  temperature  responses  of  un-heat  acclimated  5th,  50th  and  95th  percentile  sailors 
experiencing  conditions  in  disabled  submarine  while  at  rest  in  shorts  and  T-shirts. 
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Figure  7b.  Simulated  skin  temperatures  of  resting  5th,  50th  and  95th  percentile  sailors  in  disabled 

submarine  while  at  rest  in  shorts. 


Day 


Figure  7c.  Heart  rate  responses  of  5th,  50th  and  95th  percentile  sailors  experiencing  conditions  of 
the  disabled  submarine  while  at  rest  in  shorts  and  T-shirts. 
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Figure  7d.  Sweat  rates  of  the  anthropometric  5th,  50th  and  95th  percentile  sailors  at  rest  in 

simulated  disabled  USS  Salt  Lake  City. 
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Figure  7e.  Predicted  physiological  strain  levels  of  5th,  50th  and  95th  percentile  sailors 
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RESULTS 

THE  EFFECTS  OF  RATE  OF  TEMPERATURE  RISE  ON  PHYSIOLOGY 

Depending  on  the  nature  of  the  disablement,  the  submarine’s  design,  level  of  occupancy,  the 
water  conditions  and  location,  the  rate  of  inside  Tarise  may  differ  from  the  2°F  (l.ll°C)/day  rate 
measured  in  the  disabled  USS  Salt  Lake  City  tests.  The  rate  of  temperature  rise  after 
disablement,  its  cause  and  any  upper  temperature  limit  information  or  estimates  could  influence 
the  urgency  and  planning  of  the  crew’s  rescue.  Further,  the  humidity  changes  inside  also  affect 
sailor  heat  strain  and  can  be  affected  by  furnishings,  activities  and  the  type  of  disablement. 
Figures  8a  and  8b  indicate  environmental  input  conditions  for  the  SCENARIO_J  model,  and 
Figures  8c  to  8e  presents  comparisons  of  the  physiological  responses  of  un-acclimated  sailors  of 
medium  build  (e.g.,  height:  176  cm;  weight:  81kg;  percent  body  fat  (%BF):  17%,  [Hughes  and 
Horn,  2006]) )  to  these  environmental  conditions  with  3  different  rates  of  temperature  rise  (i.e., 

1,  2  and  3°F/day  (0.56,  1.11  and  1.67°C/day))  starting  from  75.2°F  (  24°C)  and  60%  RH 
(Figure  8a&b).  The  relative  humidity  and  water  vapor  pressure  (Pvap)  levels  in  the  air  (Figure 
8b)  are  assumed  to  be  driven  by  water  vapor  releases  from  the  sailors  and  modeled  to  correspond 
to  the  test  results  from  the  USS  Salt  Lake  City. 

The  relative  humidity  increases  with  time  are  very  similar  for  the  1,  2  and  3°F/day  rates  of  Ta 
rise;  reaching  100%  in  about  6  days.  However,  vapor  pressure  continues  to  increase  though 
saturated  (100%  RH)  because  the  saturated  vapor  pressure  of  water  increases  with  temperature 
enabling  the  air  to  absorb  more  water  vapor.  The  increasing  Pvap  reduces  the  vapor  pressure 
difference  between  with  water  on  the  skin  and  the  environment,  reducing  the  evaporation  rate  per 
unit  area  of  water  on  the  skin  and  causing  the  sweat  wetted  area  to  increase.  Eventually  the  skin 
is  nearly  completely  wet,  limiting  any  further  evaporative  heat  loss  increase  and  causing  the 
body  temperature  to  rise  more  rapidly  (Figure  8c). 

The  sailors  at  the  time  of  the  disablement  are  wearing  standard  uniforms  of  long  sleeved  shirt 

2  1 

and  trousers,  with  an  intrinsic  resistance  to  dry  heat  flow  of  about  0.61  clo  (0.095 °Cm'w'  ) 
between  skin  and  outer  surface  of  the  clothing.  For  comparison,  a  winter  business  suit  is  about  1 
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clo  (0.155°Cm2w  n.  Their  Tsk  would  be  about  91°F  (33°C).  When  Tsk  increases  to  about  95°F 
(35°C),  the  Sailor  will  likely  feel  some  warm  discomfort  and  it  is  assumed  they  will  then  change 
into  shorts  and  T-shirts  or  similar  light  attire,  reducing  their  clothing  insulation  to  about  0.2  clo 
(0.03°C  nrw'1)  (Fig  8c)  .  This  clothing  reduction  will  occur  after  about  one  day  with  the 
temperature  rising  slowly  at  1°F  (0.56°C)/day,  and  after  about  2  hours  for  the  2°F  (l.ll°C)/day 
and  3°F  (1.67°C)/day  conditions,  but  when  they  choose  to  reduce  their  clothing  had  little  effect 
on  final  results. 


Figure  8a.  Air  temperatures  input  for  SCENARIO_J  simulations  of  medium-built  sailor  in  a 
disabled  submarine  with  different  rates  of  slow  continuously  increasing  Ta. 
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Figure  8b.  Relative  humidity  and  water  vapor  pressure  levels  in  the  air  resulting  from  water 
vapor  releases  from  sailors  in  a  disabled  submarine  with  slowly  increasing  Ta.  Note  Pvap 
continues  to  increase  after  RH  reaches  100%  because  of  the  rising  Ta. 
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Figure  8c.  Core  and  skin  temperature  from  SCENARIO_J  simulations  of  a  resting  medium  built 
sailor  in  a  disabled  submarine  with  Ta  rising  at  three  different  rates,  starting  from  75°F  (24°C)  & 
60%  RH.  Initially  the  sailors  are  wearing  long  sleeved  shirts  and  trousers  for  24  hours  and 
change  into  shorts  and  T-shirt  shortly  after  disablement  when  reaches  about  95°F  (35°C). 
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Figure  8d.  Heart  rate  from  SCENARIO_J  simulations  of  resting  medium  built  sailors  in  disabled 
submarine  with  Ta  rising  at  three  continuous  slow  rates,  starting  from  75°F  (24°C)  &  60%  RH 


and  initially  wearing  long  sleeved  shirt  and  trousers  but  later  changing  into  shorts  and  T-shirt 
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The  heart  circulates  blood  to  satisfy  metabolic  needs  and  also  to  cool  the  core.  The  sailor  is 
assumed  to  have  a  constant  resting  level  of  activity  but  skin  blood  flow  for  cooling  increases  in 
response  to  rising  Tc  and  Tsk  until  it  reaches  its  maximum.  The  HR  to  deliver  this  flow  is 
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assumed  to  be  maintained  at  this  level  as  long  as  the  person’s  hydration  level  remains  normal 
(Figure  8d). 

The  sweat  rates  in  Figure  8e  increase  rapidly  after  Tc  regulation  begins  to  fail  because  heat 
loss  is  insufficient.  The  hypothalamus  correspondingly  increases  sweating  and  blood  flow  to  the 
skin  to  increase  heat  loss.  But  with  the  continuously  rising  air  temperature  and  humidity  of  such 
situations,  an  energy  balance  is  now  not  possible  and  both  sweat  output  and  skin  blood  flow 
reach  their  maximums  capabilities  without  a  balance  and  Tc  stability.  These  high  sweat  and  skin 
blood  flow  rates  can  be  maintained  until  heat  illness  occurs  assuming  the  sailor’s  level  of 
hydration  is  maintained.  The  Tc  and  PSI  limit  results  of  the  Figure  8c  and  8f  graphs  are 
summarized  in  Table  2.  Observe  in  Table  2,  that  the  ambient  Ta  and  humidity  at  the  beginning  of 
Tc  regulatory  failure  are  all  about  equal  (Ta=92°F  &  100%RH)  for  Ta  rise  rates  of  1,  2  and 
3°F/day.  Similarly  when  Tc  =  103°F  and  PSI  =  10,  the  environmental  parameters  are  also  equal 
at  99°F  and  100°F  respectively  with  100%  RH.  With  dehydration  sweat  rate  and  blood  flow 
could  be  reduced  resulting  in  regulatory  failure  occurring  sooner  at  lower  air  temperatures. 


Figure  8e.  Sweat  rates  from  SCENARIOS  simulations  of  medium  built  sailors  in  disabled 

submarine  with  Ta  rising. 


Figure  8f.  The  physiological  strain  index  (PSI)  from  SCENARIO_J  simulations  of  resting 
medium  built  sailors  in  disabled  submarine  with  air  temperature  rising  at  three  continuous  slow 

rates,  starting  from  75°F  (24°C)  &  60%  RH. 
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Table  2.  Summary  of  environmental  conditions  and  elapsed  time  since  disablement  to  when  Tc 


regulation  begins  to  fail,  Tc  reaches  103°F(39.5°C)  and  physiological  strain  index  (PSI)  reaches 
10  starting  from  environmental  conditions  of  75°F(24°C)  and  60%  RH. 


Rate  of  Ta  rise 

Tc  regulation  fails 

Tc  reaches  103 °F(39.5°C) 

PSI  reaches  10 

ATa/day 

time 

Ta 

RH 

time 

Ta 

RH 

Time 

Ta 

RH 

°F/day 

°C/day 

day 

°F  (°C) 

% 

day 

°F(X!) 

% 

Day 

°F(°C) 

% 

1 

0.56 

17 

92  (33.4) 

100 

23.5 

99  (37) 

100 

25.3 

100  (38) 

100 

2 

1.11 

8.5 

92  (33.4) 

100 

11.8 

99  (37) 

100 

12.7 

100  (38.1) 

100 

3 

1.67 

5.7 

92  (33.4) 

100 

7.8 

99  (37) 

100 

8.5 

100  (38.2) 

100 

EFFECT  OF  STARTING  CONDITIONS 

If  a  planned  situation,  route  or  activity  increases  the  possibility  of  disablement  and  increasing 
interior  temperatures,  the  submarine’s  temperature  and  humidity  could  be  reduced,  lengthening 
the  time  available  for  a  successful  rescue.  For  example,  Figures  9  a  -  h  show  the  effect  on 
responses  of  the  medium-built  non-heat  acclimated  sailor  to  a  pre-disablement  lowering  of  cabin 
conditions  from  75(24°C)  to  68°F  (20°C)  with  50%RH.  The  preparation  lengthened  the  safe 
rescue  window  by  40%. 


Figure  9a.  Air  temperatures  after  disablement  starting  from  pre-disablement  conditions  of  68°F 

and  50%RH. 
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Note,  from  comparing  Figure  9a  with  Table  2,  that  the  time  to  the  beginning  of  Tc  regulatory 
failure  increases  from: 

5.7  to  8  days  or  1.3  Adays  for  a  3°F/day  Ta  rise  rate; 

8.5  to  12  days  or  3.5  Adays  for  a  2°F/day  Ta  rise  rate;  and 

17  to  24  days  or  7  Adays  for  a  l°F/day  Ta  rise  rate. 

Further,  the  Aday  change  for  initial  temperatures  different  from  75 °F  are  for  a  given  dTa/dt 
are  approximately  linear  with  the  difference  in  starting  temperatures  from  75°F  and  the  Adays  for 
a  68°F  Tstart: 

Adays  (for  Tstart)  =  (Adays)/(75-68)  *  (75-Tstart), 
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Thus  for  other  starting  temperatures  before  disablement  the  safe  period  before  thermoregulation 
begins  to  fail  can  be  estimated  by  equation  2  with  Table  2  and  Figure  10a. 

i.e.:  Adays  (for  Tstart=72°F  &  2°F/day)  =  (3.5)/(75-68))*(75-2)  =  (3.5/7)*3  =  1.5day 
or  days  until  Tc  begins  to  fail  after  an  initial  temperature  of  75°F  with 
dTa/dt=2°F/day: 

Adays  (for  Tstart=72°F  &  2°F/day)  =  8. 5+1. 5  =  lOdays. 


Figure  9b.  Absolute  humidity  (also  called  humidity  ratio)  levels  after  disablement,  starting  from 

pre-disablement  conditions  of  68°F  and  50%RH. 
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Figure  9c.  Relative  humidity  levels  after  disablement,  starting  from  pre-disablement  conditions 

of  68°F  and  50%RH. 
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Figure  9d.  Clothing  insulation  changes  assumed  after  disablement  to  reduce  thermal  discomfort 
and  Tsk  with  cabin  temperatures  starting  at  68°F  (20°C)  and  increasing  steadily  at  three  different 

rates. 


Time  (day) 

Sailors  are  initially  wearing  a  long  sleeved  shirt  and  trousers  or  similar.  When  their  Tsk  reaches 
about  95°F,  they  will  feel  some  warm  discomfort  and  change  into  lighter  attire  (shorts  and  T- 
shirt)  as  indicated  in  Figure  9d,  which  causes  Tsk  to  decrease  before  rising  again  in  Figure  9e. 


Figure  9e.  SCENARIO_J  simulation  of  core  and  skin  temperatures  for  an  medium-built  sailor  in 
a  disabled  submarine  with  air  temperature  rising  from  68°F(20°C)/50%  relative  humidity  (RH)  at 

rates  of  l°F/day,  2°F/day  and  3°F/day. 
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The  time  when  Tc  begins  to  rise  rapidly  as  sailor’s  physiological  heat  loss  mechanisms 
become  insufficient  are  indicated  in  Figures  9a  and  c.  Notice  that  the  air  temperature  and 
humidity  levels  for  Tc  regulatory  failure  are  the  same  and  independent  of  the  rate  of  temperature 
rise.  These  environmental  levels  are  also  equal  to  those  indicated  in  Table  2  and  Figure  8  for 
disablements  beginning  with  higher  a  Ta  of  75°F  (24°C)  but  then  change  to  a  cooler  68°F(20°C) 
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Ta  before  disablement  occurs  give  2,  3.5  and  7  days  longer  time  to  find  and  accomplish  the 
rescue. 


Figure  9f.  SCENARIO_J  simulation  of  heart  rate  of  a  medium  built  sailor  in  disabled  submarine 
with  Ta  rising  at  3°F/day  (1.66°C/day),  2°F/day  and  l°F/day,  starting  from  68°F  (20°C)/50%RH 
and  initially  wearing  long  sleeved  shirt  and  trousers  and  changing  into  shorts  and  T-shirt  as 

indicated  in  Figure  9d. 
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Figure  9g.  SCENARIO_J  simulation  of  physiological  strain  index  (PSI)  for  sailor  medium  built 
sailor  in  disabled  submarine  with  air  temperature  rising  at  different  constant  rates,  starting  from 
68°F  (20°C)/50%  RH  and  initially  wearing  long  sleeved  shirt  and  trousers  and  changing  into 

shorts  and  T-shirt  as  indicated  in  Figure  9d. 
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Figure  9h.  SCENARIO_J  simulation  of  sweat  rate  of  medium-built  sailor  in  disabled  submarine 
with  Ta  rising  at  3°F/day  (1.66°C/day),  2°F/day  and  l°F/day,  starting  from  68°F(20°C)/50%  RH 
and  wearing  with  long  sleeved  shirt  and  trousers  before  changing  into  shorts  and  T-shirt  as 

indicated  in  Figure  9d. 
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EFFECT  OF  HUMIDITY 

As  surrounding  air  and  surfaces  become  warmer,  the  human’s  dry  heat  loss  capabilities 
decrease  and  thermoregulation  relies  increasingly  on  evaporative  heat  loss  processes  from  skin 
and  respiratory  surfaces.  The  rate  of  evaporation  depends  largely  on  the  vapor  pressure 
difference  between  water  on  the  evaporating  surface  and  in  the  surrounding  air.  Disabled 
submarine  test  measurements  (Figure  3)  and  calculations  (Figures  4  and  5)  indicate  the  water 
content  of  air  increases  steadily  and  thus  progressively  impairs  the  sailor’s  evaporative  heat  loss 
mechanisms.  If  water  vapor  in  the  air  could  be  removed  or  controlled  by  adsorption, 
condensation  or  other  means,  the  safe  stay  time  in  the  disabled  submarines  would  be  extended. 
The  following  simulations  illustrate  the  rescue  time  and  human  response  benefits  of  humidity 
control  efforts  characterized  by:  1)  maintaining  constant  absolute  humidity  or  Tdp  temperatures 
or  2)  maintaining  constant  levels  of  RH. 

1)  Maintain  Absolute  Humidity  or  Dew  Point  Temperature  at  Pre-disabled  Fevel 

Sailor  responses  (medium  built,  non- acclimated)  and  environmental  conditions  for  the 
disabled  submarine  where  absolute  humidity  can  be  maintained  constant  at  the  initial  level  are 
displayed  in  Figures  10-a-e  for  initial  conditions  of  75°F  (24°C)  and  60%RH.  For  these 
conditions  the  absolute  humidity,  dew  point  temperature  (Tdp)  and  vapor  pressure  are  0.01 13 
lb_H20/lb_dry  air,  60°F  (15.8°C)  and  13.7  mmHg  (Torr)  respectively  and  held  at  that  level  after 
disablement. 
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Figure  10a.  Disabled  submarine  with  a  constant  dew  point  of  60°F(15.8°C)  but  Ta  rising  at 
3°F/day  (1.66°C/day),  and  2°F/day  (l.ll°C/day),  starting  from  75°F  (24°C)/60%  RH. 
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At  the  time  of  disablement  sailors  are  wearing  long  sleeved  shirt  and  trousers  but  change  into  T- 
shirt  and  shorts  as  indicated  by  rapid  drop  in  Tsk  (Fig  10c). 


Figure  10b.  Relative  humidity  and  vapor  pressure  (Pvap)  of  water  in  disabled  submarine  with  a 
constant  dew  point  of  60°F  (15.8°C)  but  Ta  rising  at  3°F/day  (1.66°C/day),  and  2°F/day 
(1.1  l°C/day),  starting  from  75°F  (24°C)/60%  RH. 
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Figure  10c.  SCENARIO_J  simulation  of  core  and  skin  temperature  of  medium  built  sailor  resting 
in  a  disabled  submarine  with  constant  dew  point  of  60°F  (15.8°C)  but  with  Ta  rising  at  3°F/day 
and  2°F/day  (1.66°C/day  and  l.ll°C/day),  starting  from  75°F  (24°C)/60%  RH. 


40 

39 

38 

37 

36 

35 

34 

33 

32 

31 

30 


n 

o 

CD 

0) 

3 

o. 

H 

CD 

3 

TJ 

CD 

0) 

c 

CD 


i/> 

?r 


Time  (Day) 


In  comparison  to  uncontrolled  cabin  humidity,  the  safe  stay  times  with  humidity  held 
constant  at  initial  levels  (Figure  10c)  are  much  longer  than  those  with  rising  humidity  (Figure  8c 
and  Table  2).  The  constant  low  water  vapor  pressure  in  the  air  enables  sweat  to  evaporate  and  be 
an  effective  cooling  mechanism  longer  and  compensates  the  higher  cabin  temperatures. 

Figure  lOd.  Simulation  of  heart  rate  (dotted  line)  and  physiological  strain  index  (PSI,  solid  lines) 
of  a  medium-built  sailor  resting  in  a  disabled  submarine  with  constant  dew  point  of  60°F 
(15.8°C)  but  Ta  rising  at  3°F/day  and  2°F/day  (1.66°C/day  and  l.ll°C/day),  starting  from 

75°F(24°C)/60%RH. 
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Figure  lOe.  Simulation  of  sweat  rate  responses  of  medium  built  50th  percentile  sailor  in  disabled 
submarine  with  constant  dew  point  of  60°F  (15.8°C)  but  Ta  rising  at  3°F/day  and  2°F/day, 
starting  from  75°F  (24°C)/60%  RH  and  wearing  long  sleeved  shirt  and  trousers  but  changed  into 

shorts  and  T  shirt  quickly  after  disablement. 
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2)  Constant  Relative  Humidity 

The  effects  on  un-acclimated  medium  built  sailors  to  constant  RH  but  rising  Ta  are  illustrated 
in  the  following  graphs  (Figures  lla-f).  As  in  previous  simulations  the  occupants  are  resting  and 
initially  wearing  a  long  sleeved  shirt  and  trousers  but  change  into  shorts  and  T-shirt  soon  after 
disablement.  All  Figure  1 1  simulations  are  for  Ta  rising  at  the  rate  of  l°F/day  (0.56°C/day). 

Figure  11a.  Simulation  of  environment  in  disabled  submarine  with  air  temperature  (Ta)  rising  at 
l°F/day  (0.55°C/day)  from  an  initial  75°F  (24°C)  and  60%  relative  humidity  (RH)  condition 
with  RH  held  constant  after  a  rise  from  the  starting  level. 
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Figure  1  lb.  SCENARIO_J  simulation  of  core  temperature  of  a  resting  medium  built  sailor  in  a 
disabled  submarine  with  relative  humidity  (RH)  held  constant  after  rise  from  starting  level  and 
with  air  temperature  (Ta)  increasing  at  l°F/day  (0.55°C/day)  from  initial  75°F  (24°C)  and  60% 
RH  condition.  After  dayl,  the  sailors  are  wearing  shorts  and  T-shirt. 
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Sailors  are  wearing  long  sleeved  shirt  and  trousers  at  disablement  but  change  into  shorts  and  T- 
shirt  as  indicated  by  rapid  drop  in  Tsk  (Figure  1  lc). 


Figure  lie.  SCENARIO_J  simulation  of  skin  temperature  of  a  resting  medium  built  sailor  in  a 
disabled  submarine  with  relative  humidity  (RH)  held  constant  at  indicated  levels  and  with  Ta 
rising  at  l°F/day  (0.55°C/day)  from  initial  75°F  (24°C)  and  60%  RH  condition. 


Time  (day) 
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Figure  1  Id.  Heart  rate  response  of  resting  50th  percentile  sailor  in  a  disabled  submarine  with  RH 
held  constant  and  with  Ta  rising  at  l°F/day  (0.55°C/day)  from  initial  75°F  (24°C)  and  60%  RH 

condition. 
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Figure  1  le.  Simulation  of  PSI  responses  of  50th  percentile  sailor  resting  in  a  disabled  submarine 
with  RH  held  constant  after  rising  from  starting  level  and  with  Ta  rising  at  l°F/day  (0.55°C/day) 

from  initial  75 °F  (24°C)  and  60%  RH  condition. 
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Figure  1  If.  Simulation  of  sweat  rates  of  resting  50th  percentile  sailor  in  disabled  submarine  with 
RH  held  constant  after  rising  from  starting  level  and  with  Ta  rising  at  l°F/day  (0.55°C/day)  from 

initial  75°F  (24°C)  and  60%  RH  condition. 
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The  evironmental  and  time  limit  values  of  Figure  lib  are  summarized  in  Table  3.  The  table 
gives  the  results  for  a  temperature  rise  rate  of  l°F/day.  For  other  slow  rates  of  temperature 
change  (X  °F/day)  the  elapsed  time  (Atx)  of  an  event  can  be  estimated  by: 

Atx  =  (Atx  for  l°F/day)/(X°F/day). 

Table  3.  Summary  of  values  from  Figure  1  la,b  and  e  for  constant  relative  humidity  (RH)  and 
time  where  core  temperature  (Tc)  regulation  begins  to  fail,  Tc  =  38.5°C,  Tc  =  39.5°C,  and  PSI  = 
10  with  submarine’s  interior  air  temperature  rising  at  l°F/day  starting  from  75°F  (24°C)  and 
60%  RH  and  sailors  resting  in  shorts  and  T-shirt. 


Tc  regulation  begins 

to  fail 

Tc  reaches 

10TF 

(38.5°C) 

Tc  reaches 

103°  F 

(39.5°C) 

Tc  reaches 

105°  F 

(40.5°C) 

PSI 

reaches  10 

RH 

time 

Ta 

Tc 

time 

Ta 

time 

Ta 

time 

Ta 

time 

Ta 

(%) 

(day) 

°F(°C) 

°C 

(day) 

°F(°C) 

(day) 

°F(°C) 

(day) 

°F(°C) 

(day) 

°F(°C) 

100 

17.5 

92(33.4) 

37.08 

21.5 

96.8(37) 

23.3 

98.6(37) 

25.25 

100.4(38) 

25.4 

100(38.1) 

90 

20.2 

95.4(35.2) 

37.12 

24.1 

99(37.2) 

26 

101(38.4) 

27.8 

103(39.5) 

28.1 

103(39.5) 

80 

23.5 

98.6(37) 

37.22 

27 

102.2(39) 

29 

104.2(40.1) 

31 

106.2(41.2) 

31.1 

106.2(41.2) 

70 

26.8 

102(38.9) 

37.22 

30.4 

105.4(40.8) 

32.4 

109.9(43.2) 

34.3 

109.8(43.2) 

34.4 

109.8(43.2) 

60 

30.6 

105.6(40.9) 

37.25 

34.1 

109.4(43) 

36.2 

111.6(44) 

38.25 

113.5(43.3) 

38.4 

113.5(45.4) 

COMPARISON  OF  AIR  TEMPERATURE  LIMITS  AND  THE  TOLERANCE  TIME 
AFTER  DISABLEMENT  AS  PREDICTED  BY  SCENARIOS  AND  SUB_ICDA 
MODELS 

The  predicted  Ta  and  times  to  reach  the  beginning  of  Tc  regulation  failure  and  also  for  Tc  to 
reach  101°F  (38.5°C)  and  103°F  (39.5°C)  are  displayed  in  Figure  12  for  the  SCENARIO_J  and 
Sub-ICDA  models.  The  predictions  are  for  3  rates  of  temperature  rise  and  2  starting  points.  The 
air  temperatures  at  which  Tc  fails,  reaches  101  °F  (38.5°C)  and  103°F  (39.5°C)  are  independent 
of  the  Ta  rise  rates  and  starting  conditions  for  both  models.  This  is  because  the  three  rates  (1,2 
and  3°F/day)  of  temperature  rise  are  very  slow,  enabling  the  sailor’s  physiology  and  body 
temperatures  to  adjust  and  be  in  a  quasi-steady  condition  with  the  environment. 

The  Ta  levels  for  the  Tc  limits  indicated  in  Figure  12  are  for  disabled  submariners  with 
uncontrolled  humidity  that  rises  and  reaches  100%  RH  before  the  sailor’s  Tc  begins  to  rise 
rapidly  from  progressive  insufficient  heat  loss  and  Tc  regulation  failure.  The  Ta  at  the  Tc  limits 
by  the  two  models  are  very  similar.  The  Sub_ICDA  model  consistently  predicting  the  limits  to 
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occur  about  0.6°C  (1°F)  lowerTa  ,  except  at  beginning  of  Tc  regulation  failure  that  is  about  1.6°C 
(3°F)  lower  than  the  Ta  SCENARIO_J  model  predicts. 


Figure  12.  Comparison  of  air  temperatures  (Ta)  where  core  temperature  (Tc)  regulation  begins  to 
fail,  Tc  reaches  38.5°C  and  39.5°C  (101. 3°F  and  103. 1°F)  after  dissablement  as  predicted  by 
SCENARIO_J  (upper  graph)  and  Sub_ICDA  (lower  graph)  models. 
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A  comparison  of  PSI  levels  predicted  by  the  models  at  the  Tc  limits  and  times  are  displayed 
in  Figure  13.  The  agreement  is  again  consistent  with  Sub  _ICDA  predicting  about  1.7  PSI  units 
less  strain  at  the  38.5°C  and  39.5°C  Tc  limits  and  0.5  units  less  at  the  beginning  of  Tc  failure. 


Figure  13.  Comparison  of  physiological  strain  (PSI)  predictions  by  SCENARIO_J  (upper  graph) 
and  Sub_ICDA  (lower  graph)  at  equal  environmental  conditions. 
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The  Ta ,  time  and  PSI  values  of  Figures  12  and  13  for  the  two  modeling  systems  are 
summarized  in  Table  4.  As  noted  above,  the  Ta’s  and  their  occurrence  times  predicted  by  the 
models  are  very  similar,  particularly  for  when  Tc  reaches  101°F  (38.5°C)  and  103°F  (39.5°C). 
The  Ta  are  generally  a  little  lower  and  occur  sooner  for  the  SubJCDA  response  predictions. 
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The  coincident  PSI  levels  are  also  lower  for  Sub_ICDA  at  the  Tc  limits.  The  SCENARIO_J 
model  has  the  more  complete  and  developed  thermal  physiological  control  mechanisms,  and  may 
be  the  more  accurate  predictor  for  this  situation.  But  the  differences  are  small  and  the  response 
patterns  are  similar  making  both,  or  their  average,  good  guidelines  for  rescue  planning. 


Table  4.  Summary  of  time  elapsed  since  disablement,  cabin  Ta,  and  sailor’s  PSI  values  where:  Tc 
regulation  begins  to  fail,  TC=101.3°F  (38.5°C)  and  TC=103°F  (39.5°C)  as  predicted  by 
SCENARIO_J  and  Sub_ICDA  models  for  resting  sailors  in  shorts  and  T-shirts.  The  disabled 
submarine’s  humidity  100%  RH  before  Tc  regulation  begins  to  fail  and  remains  at  100%  RH 

thereafter . 


Tc  regulation  begins  to  fail 


Tc=101.3°F(38.5°C) 


model 

dTa/dt 

starting  at 

time 

Ta 

Ta 

PSI 

time 

Ta 

Ta 

PSI 

time 

Ta 

Ta 

PSI 

°F/day 

°FCC)  /%RH 

day 

°F 

°C 

day 

°F 

°C 

day 

°F 

°C 

ScenarioJ 

1 

75(24)/60 

17.5 

92.1 

33.4 

0.9 

21.5 

96.6 

35.9 

6 

23.4 

98.4 

36.9 

8 

2 

8.4 

92.1 

33.4 

0.9 

10.8 

96.8 

36 

5.8 

11.8 

98.8 

37.1 

7.9 

3 

5.8 

92.67 

33.7 

1 

7.25 

97 

36.1 

6 

7.8 

98.8 

37.1 

7.8 

1 

68(20)  /50 

24 

92.1 

33.4 

0.8 

28.5 

96.8 

36 

5.9 

30.4 

98.4 
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2 

12 

91.9 

33.3 
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97 
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6 
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99 
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3 

8 
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33.3 

0.9 

9.75 

97.3 

36.3 

6 

10.3 

98.8 

37.1 

8.9 

96.9 


36.1 


Sub-lcda  1  75(24)/60  13.5  88.9  31.6  0.4 

20 

95.3 

35.2 

4.3 

22  97.3  36.3 

6.3 

2  6.9  88.7  31.5  0.45 

10.2 

95.7 

35.4 

4.3 

11  97.6  36.45 

6.4 

3  4.8  89.4  31.9  0.5 

6.75 

95.5 

35.3 

4.2 

7.4  97.5  36.4 

6.3 

1  68(20) /50  20.5  88.7  31.5  0.4 

27.2 

95.5 

35.3 

4.4 

29.2  97.5  36.4 

6.4 

2  10.4  88.9  31.6  0.5 

13.7 

95.4 

35.2 

4.3 

14.7  97.5  36.4 

6.2 

3  7  89.1  31.7  0.4 

9.2 

95.7 

35.4 

4.3 

9.8  97.5  36.4 

6.5 

95.5 


35.3 


Escape 

Previous  model  simulations  are  based  on  conditions  for  the  resting  sailor  wearing  shorts 
and  T-shirt  or  similar  clothing  to  survive  longer  periods  of  time  in  a  disabled  submarine. 
However,  escape  from  the  disabled  submarine  may  necessitate  wearing  an  escape  suit. 
Unfortunately  for  these  warm  conditions,  the  impermeable  escape  suit  eliminates  possible 
evaporative  cooling  from  sweating;  causing  the  sailor’s  Tc  to  rise  more  rapidly  and  shortens  the 
time  available  for  safe  crew  rescue.  Physiological  responses  to  the  disabled  thermal  conditions 
were  estimated  by  the  SCENARIO_J  model  for  the  50th  percentile  sailor  donning  the  escape  suit 
when  the  Ta  reaches  86°F  (30°C)  and  also  when  thermal  regulation  begins  to  fail  at  an  air 
temperature  of  92°F  (33°C)  with  100%  RH. 

The  Tc  responses  are  displayed  in  Figure  14  for  1,  2  and  3  °F/day  rates  of  Ta  rise.  Donning 
the  escape  suit  at  92°F  conditions  causes  Tc  of  the  resting  sailor  to  reach  the  unsafe  103°F 
(39.5°C)  level  in  about  2  hours.  If  the  suit  is  put  on  earlier  when  the  interior  temperature  is  86°F 
(30°C),  the  time  available  before  reaching  103°F  (39.5°C)  increases  to  about  1  day  irrespective 
of  the  rate  of  temperature  rise.  Figure  14  strongly  indicates  that  the  wearing  of  the  escape  suit 
should  be  delayed  until  the  escape  begins,  and  avoided  if  possible  in  such  thermally  challenging 
situations. 
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Figure  14.  Core  temperature  responses  to  donning  an  escape  suit  at:  1)  Ta  =  86°F  (30°C),  2)  Tc 
regulation  begins  to  fail,  and  3)  not  donning  the  escape  suit.  Cabin  temperature  is  rising  at 
l°F/day  (upper  graph),  2°F/day  (middle),  and  3°F/day  (lower  graph),  starting  from  75°F  (24°C) 
and  60%  RH.  Humidity  is  uncontrolled  and  reaches  100%  RH  as  indicated  at  top  of  each  graph. 

Time  axes  of  graphs  have  different  day  lengths. 
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DISCUSSION 


The  sailor  response  predictions  were  made  with  reliable  USARIEM  thermo-physiological 
models  developed  and  validated  for  warm  environments.  The  simulations  estimate  human 
responses  to  a  loss  of  power  in  an  intact  vessel  resulting  in  a  slow  steady  increase  of  air 
temperature  and  humidity  in  the  submarine’s  occupied  compartments.  In  some  disablements 
equipment  damage,  water  leaks  and  personal  injuries  may  impose  additional  stress  and  strains  on 
the  sailors  altering  the  assumed  modeled  conditions  and  resting  activities  while  waiting  for  the 
rescue.  Increasing  the  average  metabolic  activity  of  the  crew  will  in  general  shorten  the  time 
available  for  safe  successful  rescue.  Further,  the  study  assumes  the  submarine’s  occupied 
compartment  temperatures  can  increase  without  limit.  However,  depending  on  the  situation,  air 
temperature  upper  limits  will  likely  develop  limiting  potential  thermal  stress.  Because  of  these 
possible  unknowns  in  modeling  the  responses  of  the  crew  and  also  the  condition  of  the 
submarine  after  disablement,  the  results  of  this  study  should  be  considered  as  rational  guidelines 
and  be  compared  to  observations  and  measurements  when  possible  as  they  become  available  or 
are  discovered  from  applicable  sea  and  laboratory  experiences. 

The  sailor’s  responses  to  the  slow  rise  of  temperature  in  the  disabled  submarine  have  been 
simulated  for  a  wide  range  of  conditions.  The  results  have  similar  repeatable  patterns.  The 
gradual  temperature  and  humidity  increase  permits  the  sailors  thermoregulatory  system  to  be  in  a 
quasi  equilibrium  state  with  the  environment.  This  is  fortunate  and  enables  the  environmental 
parameters  important  for  heat  exchange  to  be  reasonable  and  reliable  indicators  of  the  sailor’s 
thermo-physiological  state.  Figure  15,  developed  from  simulation  results,  based  on  Tables  2,  3 
and  4  summarizes  the  relationship  between  the  environmental  temperature  and  humidity  with 
quasi-steady  states  of  body  Tcfor  the  resting  sailor  wearing  shorts  and  T-shirt  or  similar  clothing. 
The  body  temperatures  are  independent  of  the  rate  of  temperature  for  rates  <6°F/day.  Thus  the 
easily  measured  air  temperature  and  humidity  of  the  occupied  spaces  is  a  good  indicator  of  the 
thermal  health  state  of  the  occupants.  And  they  can  both  be  readily  and  reliably  measured.  The 
time  for  the  average  sailor  to  reach  the  beginning  of  Tc  regulation  failure,  and  or  Tc  =  101°F 
(38.5°C),  Tc  =  103°F  (39.5°C)  and  Tc  =  105°F  (40.5°C)  can  be  estimated  from  initial 
temperature,  the  rate  of  Ta  rise,  and  the  Ta  for  the  event  from  Figure  15  and  equation  [4],  as 
explained  in  example  1. 
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Figure  15.  The  air  temperature  and  humidity  levels  for  indicated  body  core  temperatures  of 
resting  crew  in  a  disabled  submarine  with  slowly  increasing  temperature. 


Humidity  effects  on  cabin  temperature  limits  for  resting 
sailors  in  shorts  and  T-shirt 
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Elapsed  time  since  disablement  for  Tc  to  reach  a  level  indicated  in  Figure  5 
Elapsed  time  =  (Ta@Tc  -  Ta@disablement)/( average  Ta  rise) 


Example  1.,  Time  to  when  Tc  regulation  begins  to  fail  with  RH=75%  and  an  average  Ta  rise  of 
2.2°F/day,  starting  from  disablement  cabin  air  temperature  of  71°F.  From  figure  15  at  Tc  fails  & 
75%RH,  Ta=101.5°F.  Then,  substituting  into  equation  4,  the  elapsed  time  =  (101.5- 
71)/2.2=13.86  days  and  Tc  regulation  failure  begins  at  37.22°C  (99°F)  (Figure  15). 

Example  2.  When  Ta  is  98°F  and  RH  is  95%,  then  quasi-steady  state  of  Tc  is  approximately 
101.3  °F  (Figure  15). 
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Air  Temperature  Ta  (°C) 


CONCLUSIONS 


Thermo-physiological  responses  of  crew  members  to  the  gradual  but  unrelenting 
accumulation  of  heat  and  humidity  possibilities  in  a  disabled  submarine  were  simulated  with 
USARIEM  computer  modeling  tools.  The  modeling  indicates  that  the  resting  medium  body  size 
(50th  percentile)  sedentary  and  lightly  clad  sailor,  waiting  for  rescue,  would  begin  to  experience  a 
rapid  rise  in  body  temperature  and  thermo-regulatory  failure  at  an  air  temperature  (Ta)  about 
92°F  (33°C)  together  with  humidity  near  saturation  (95-100%  RH).  The  humidity  rise  was 
assumed  to  be  driven  by  water  vapor  evaporating  from  the  crew.  From  that  point  on,  it  will  take 
6  more  days  with  a  1°F  (0.6°C)/day  rate  of  Tarise  for  the  individual’s  Tc  to  reach  the  unsafe 
103°F  (39.5°C)  level,  and  3  and  2  days  with  2°F  (1.1°C)  and  3°F  (1.7°C)/day  rise  conditions, 
respectively.  Heat  acclimation  status  appears  to  have  minimal  effect  on  final  body  temperatures 
and  rescue  time  limits.  Short-  lean  sailors  will  be  slower  in  reaching  unhealthy  body 
temperatures  and  sailors  with  large  body  size  will  arrive  at  these  levels  sooner.  The  time 
available  for  the  rescue  can  be  lengthened  by  cooler  and  dryer  initial  pre-disablement  conditions. 
Also  reducing  or  slowing  the  rise  of  humidity  has  a  strong  beneficial  effect  on  thermoregulation 
with  rising  warm  temperatures. 

At  92°F  (33°C)  with  high  humidity  the  individual  is  relying  on  sweat  evaporation  to  carry 
away  body  heat.  Donning  an  impermeable  escape  suit  eliminates  the  evaporative  cooling  from 
sweating,  causing  Tc  to  rise  rapidly.  At  92°F  (33°C)  conditions  in  the  escape  suit,  Tc  reaches 
103°F  (39.5°C)  in  about  2  hours.  If  the  suit  is  put  on  earlier  when  the  interior  temperature  is  86°F 
(30°C),  the  time  available  before  reaching  103°F  (39.5°C)  increases  to  about  1  day  irrespective 
of  the  rate  of  temperature  rise.  This  strongly  indicates  that  the  wearing  of  the  escape  suit  should 
be  delayed  until  the  escape  begins,  and  avoided  if  possible  in  such  situations. 


RECOMMENDATIONS 

The  project  developed  a  useful  graph  (Figure  15)  indicating  cabin  temperature  and  RH 
conditions  that  correspond  to:  1)  beginning  of  thermoregulatory  failure  and  rapid  rise  of  Tc,  2) 
Tc  =  101. 4°F  (38.5°C),  3)  Tc  =  103°F  (39.5°C),  and  4)  Tc  =  105°F  (40.5°C).  As  Ta  and  RH  can 
be  readily  measured  by  the  crew,  the  graph  may  help  estimate  the  thermal  health  status  of  the 
sailors  and  guide  rescue  activities  after  disablement.  Further,  the  graph  and  accompanying  Tables 
2,  3  and  4  can  be  used  to  estimate  the  expected  time  to  reach  these  various  health  states. 
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APPENDIX  A 


Listing  of  Sub_ICDA  source  code 

/* 

Origin:  BBMD  @  USARIEM 
Authors:  Berglund  and  Yokota 
Date:  3/31/08 

File:  Sub_ICDA_HRk_360_3_31_08 . c 

This  simulates  the  thermal  physiological  responses (Tc, Tsk, sweat,  skin 
diffusion  and  repiratory  waterLoss , skin  wettedess , HR,  PSI)  from  inputs  of 
metabolism,  environmental  parameters (Ta, RH, V,  Air  Volume) and  constant 
personal  properties (ht, wt  %Bfat)  and  clothing. 

The  constant  personal  properties , and  initial  physiological  values  are  entered 
through  the  keyboard.  The  metabolism  and  other  parameters  are  entered  at 
different  simulation  times  are  read  from  a  textfile.  The  interval  (in 
seconds)  between  measurements  is  adjustable.  Enter  0  for  the  default  interval 
of  300sec (5min) . 

The  measured  values  (min, met)  are  read  from  a  txt  file.  See  Subln2.txt  for  an 
example  of  input's  form  and  spacing (consecutive  10  digit  spaces  for  each  min 
of  measurements) . 


The  output  is  displayed  on  screen  and  written (appended) to  output.txt  file. 
The  output  of  predicted  responses  occur  after  each  integration  step. 

The  program's  integration  step  is  adjustable.  Enter  0  for  the  default 
integration  step  interval  of  60sec.The  default  value  is  a  good  choice,  for 
slow  changing  situations  interval  step  can  be  longer  and  a  shorter  interval 
for  responses  to  rapidly  changing  stituations. 

Metabolism  is  expressed  by  dimenionless  relative  metabolism  term (met)  where 
met=actual  metabolism/resting  metabolism. 


intrinsic  clo  is  thermal  resistance  from  skin  to  outer  layer  of  clothing  in 
clo  units. 

MRT=  mean  radiant  temperature 
V=air  speed  m/s. 

*/ 


#include  <assert.h> 
#include  <stdio.h> 
#include  <stdlib.h> 
#include  <math.h> 
#include  <ctype.h> 


#def ine 
#def ine 
#def ine 
#def ine 
#def ine 
#def ine 
#def ine 
#def ine 
#def ine 


Pbar  760. 

//atmospheric  pressure  mmHg 

TTSK  33.7 

// 

"C 

TTCR  36.8 

// 

°C 

TTBM  36.49 

// 

TTBM= . 9*TTCR+ . 1*TTSK  =36.49 

CSW  170. 

// 

g/ (h  mA2  °C) 

CDIL  50. 

// 

CDIL  200  super  athlete,  50  average  person 

CSTR  .5 

// 

1/°C 

SKBFN  6.3 

// 

liters/ (h  mA2 ) 

Skbfmax 

90  . 

//  conservative  could  be  higher  for  fit  person 
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#def ine 
#def ine 
#def ine 
#def ine 
core 
#def ine 
muscle 
#def ine 
#def ine 
#def ine 
#def ine 
#def ine 


Skbfmin  2 . 
CMIN  5.28 
Clave  3200. 
F02ExtCr  .10 

F02ExtMu  .80 

BL02cap  .21 
02equiv  2.93 
CB  1.163 
WpMet  58.2 
LR  2 .2 


//Liters/ (h  mA2) 

//  w/ (mA2  °C) 

//  mlBlood/ (mA2*min) 

//fraction  of  02  extracted  from  blood  flowing  through 

//fraction  of  02  extracted  from  blood  flowing  through 

/ /ml02/mlBlood 
//ml02/ (w*min) 

//  wh/ (L  °C) 

//  w / (m2*met) 

//Lewis  relation  °C/mmhg 


void  OutputHeader ( ) ; 


double  SatVapPres (double  T) ;  //  in  Torr  =  mmhg 

double  Clot(int  code,  double  V);  //total  clo  from  skin  to  environment 
double  VRes(int  code,  double  V);  //total  vapor  resistance  from  skin  to 
environment 

double  Convection (double  met,  double  V) ; //convective  heat  transfer  coef 
double  Respiration (double  Rm,  double  Ta,  double  Pa, double*  Wres); 
//respiratory  heat  loss 

double  SkinBloodFlow (double  Tc,  double  Tsk); 
double  Shiver (double  Tc,  double  Tsk, double*  Bfat); 

double  DuBoisSkinArea (double  Wt,  double  Ht) ;  //surface  area  of  body 

double  02pkg (double  Bf,  double  HR); 

double  BflowCore (double  mtot, double  HRRk); 

double  BflowMuscle (double  mtot, double  Bfc); 

double  StrokeVolume (double  Adu,  double  HRrest, double  Tsk); 

double  BreathRate (double  met, double  Wt, double  Adu)  ; 

double  MetT (double  HR, double  HRrest, double  T) ; 

double  HeartRate (double  met, double  T,  double  HRrest);  //from  ICDA 

double  Alpha (double  Skbf ) ;  //  fraction  of  body  that  is  ski 

double  Sweat (double  Tc,  double  Tsk, double  alpha); 

double  Core (double  Tc,  double  Tsk, double  met, double  Res, double*  hfcsk,  double 
skbf, double*  Rm,  double*  Bfat); 

double  Skin (double  wet, double  Dry,  double  Emax,  double*  hfcsk); 

double  MoranPSI (double  Tc, double  Tco, double  HR, double  HRrest);  //heat  stress 
index 


int  main  (void) 

{ 

int  step,  code; 

double  Age, Ta,Tao,Ta  rise, MRT, RH, RHo, Pa, Psa, Psk, Tel, To, V, Vps , ma, HRkmax; 
double  HRRk, he, he, hr, met, Metnew, Ht , Wt , Adu, Wres , drip,  DT,  DTs ,  DTM,  DTMP,  Res ; 
double  DTMH, ExTime, MeasTime; ReadTime, Print Time, rhoa, clos , cloDisSub, Emax; 
double  CEvapWLL, V02max; wet, Tdpa, Tsk, Tsko, Tc, Tco, Rm, clo, FACL, Dry, PSIc, BMI 
double  Esk,Edif,  Skbf , HRrest, HR, alpha, Regsw, Ersw, Wi, Wo, Wos, Wdif, Wsk  evap 
double  Wevap, heatFlowCoreToSkin, HSCR, HSSK, TCCR, TCSK, Rcl, Rbound, Rclt,  Rpcl 
double  Rpbound, Rpclt, Bf , Bfc, Bfm, Bfat, SV, CEvapWL, CdripL  time, TIM, MTIM; 
double  Duration, CumulativeWaterLossPsqm, CWLL;  CumulativeWaterLoss; 
doubleRateOf Total WaterLossPsqm, RateOf TotalWaterLoss  ; 
char  inName[64],  outName[64]; 

FILE  *pif ile=NULL;  //file  pointer  to  inputfile 

FILE  *pof ile=NULL;  //file  pointer  to  outputfile 

char  *inputf ilename=inName; 
char  *outputf ilename=outName; 
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//  get  file  names  from  user 

printf ( "Enter  name  of  inputfile:  xxxx.txt  "); 
scanf  ("%63s",  inName) ; 

printf ( "Enter  name  of  outputfile (<30char) :  xxxxx.txt  "); 
scanf  ("%63s",  outName) ; 

TIM=0 ;  MTIM=0;  step=0;  time  =0;  HR=0;  HRrest=0;  DTs=0;  ExTime=0; 
ReadTime=0;  MeasTime=0;  PrintTime=0;  Wo=0;  Wos=0;  Wi=0;  drip=0; 

Wo=0;  Wos=0;  Wi=0;  CumulativeWaterLossPsqm=0 ;  CumulativeWaterLoss=0 ; 
CWLL=0 ;  CEvapWL=0 ;  PSIc=0;  Regsw=0;  CEvapWL=0;  CdripL=0; 

printf ( "Enter  time-step  for  discrete  integration  in  seconds (0  for  default 
of  60  seconds) :DTs  "); 
scanf ("%lf",  &DTs )  ; 
if  (DTs<=0 ) 

DTs=60;  //  default  time  step  in  seconds 

DT=DTs/ (60*60) ;  //time  step  DT  in  hrs 

printf ( "Enter  interval  between  measurements  in  minutes ( interval>=DTs ) (0 
for  default  of  (60min) :  "); 

scanf ("%lf ", &DTM) ; 
if  ( DTM<=0 ) 

DTM=60;  /*time(min)  between  measurements  default  value  of  60  min 
between  measured  input  */ 

DTMH=DTM/ (60);  //time  between  measurements  DTMH  in  hrs 

printf ( "Enter  starting  intrinsic  clo (sailor  fatigues=.6)  :  clo  "); 

scanf ("%lf", &clos) ; 

code=0 ; 

FACL=1+ . 2*clos;  //estimate  of  outside  area  of  clothed  person 
Rcl= . 155*clos; 

Rpcl=  . 153181*clos;  /*  Berglund  1985, for  icl=Rcl/ (LR*Rpcl) = . 45; 

Rpcl= . 155*clo/ (LR*icl) =. 155*clo/ (2. 2*. 45)  */ 

clo=clos ; 

printf ( "Enter  DisSub  intrinsic  clo ( shorts , Tshirt, sandals= . 2  )  :  clo  "); 

scanf ("%lf", ScloDisSub) ; 

printf ( "Enter  subject's  weight ( kg) , height (m)  (or  0's  if  unknown):  Wt  Ht 
");/*  inserts  default  values  of  82kg, 1.77m  */ 
scanf ("%lf%lf", &Wt, &Ht) ; 
if  (Wt<=0) 

Wt=82 ; 
if  (Ht<=0) 

Ht=l .77; 

printf ( "Enter  subject's  %BFat  (or  0's  if  unknown):  Bfat  ");/*  inserts 

default  values  of  15%  */ 

scanf ("%lf", &Bfat) ; 
if  (Bfat<=0) 

Bfat=15; 

printf ( "Wt=%7 . 3f  kg,  Ht=%7.3f  m,  Bfat=%8.4f  \n",  Wt,Ht,Bfat); 
Adu=DuBoisSkinArea (Wt, Ht) ;  / /mA2 


printf ( "Enter  subject's  resting  HR, age  and  V02max (ml/min/kg) (or  0's  if 
unknown) :HRrest  Age  V02max  ");/*inserts  default  values  of  70,23,51.4  for 
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zero ' s  * / 

scanf ( "%lf %lf %lf " , &HRrest, &Age, &V02max) ; 
if  (HRrest<=0) 

HRrest=7  0 ; 
if  (Age<=0) 

Age=23; 

if  (VO2max<=0)  //ml/ (min  kg) 

V02max=51 . 4 ; 

printf ( "Enter  initial  physiology  (or  0's  if  unknown) :  Tsk  Tc  HR  met  ");/* 

inserts  default  values  if  Tsk, Ter  are  zero's  */ 
scanf ( "%lf%lf %lf %lf " , &Tsk, &Tc, &HR, &met) ; 
if  (Tsk  <=0 ) { 

Tsk=33 ; 

Tsko=Tsk; 

} 

else 

Tsko=Tsk; 

if  (Tc  <=0) { 

Tc=36 . 9; 

Tco=Tc; 

} 

else 

Tco=Tc; 

if  (HR<=0 ) 

HR=HRrest; 
if  (met<=0) 
met=l ; 

Rm=  met*58.2; 

HRkmax=2  2  0 -Age ; 

HRRk=HR/ HRkmax ; 

BMI=Wt/ (Ht*Ht) ; 
wet= .06; 

printf ( "Enter  initial  environment  (or  0's  if  unknown) :  Ta  RH  V  ");/* 

inserts  default  values  if  Ta,RH  V  are  zero's  */ 
scanf ( "%lf %lf %lf " , &Ta, &RH, &V) ; 
if  (Ta  <=0 ) { 

Ta=24 ; 

Tao=Ta; 

} 

else 

Tao=Ta; 

if  (RH  <=0 ) { 

RH=60 . ; 

RHo=RH; 

} 

else 

RHo=RH; 

if  (V<=0) 

V=. 05; 

MRT=Ta; 
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printf ( "Enter  air  volume  per  sailor  and  initial  rate  of  Ta  rise  (or  0's 
if  unknown) :  Vps  Ta^rise  ");/*inserts  default  values  if  Vps,Ta  rise  are 

zero ' s  *  / 

scanf ( "%lf%lf" , &Vps , &Ta_rise) ;  //units  mA3,  C/day 

if  (Vps  <=0) 

Vps=360;  //mA3 

if  (Ta_rise<=0) 

Ta_rise=l . 1 1 ;  //  2F/day  =  1 . llC/day=l . 11/24  C/h 

Pa=RH*SatVapPres (Ta) /100;  //  ambient  vapor  press  Torr 

Tdpa=4030 .183/ (18.6686-log(Pa))-235; 

Wi= . 622*Pa/ (Pbar-Pa)  ;  //  initial  humidity  ratio  kgw/kga 

Wo=Wi ; 

rhoa= (Pbar-Pa) *133/ ( (Ta+273) *8315/28 . 97) ;  //density  of  dry  air  kga/mA3 
ma=rhoa*Vps;  //kg  of  air 

printf ( "rhoa=%7 . 3f  kg/mA3,  ma=%7.3f  kga,  Wi=%8.4f  \n",  rhoa,ma,Wi); 

printf ( "Enter  duration (hr) (0  for  default=16days=384h) :  ");  /*  may  wish 

to  delete  duration  for  real  applications  */ 
scanf ( "%lf" , SDuration) ;  //  duration  in  hours,  hr 

if  (Duration<=0 ) 

Duration=384 ;  //16days 

printf ( "Enter  interval  in  min  between  printouts (0  for  default  of 
(360min=6hr) :  "); 

scanf ("%lf", &DTMP) ;  / /DTMP  in  min 

if  ( DTMP<=0 ) 

DTMP=360;  / / 6  hours 

pifile=fopen (inputf ilename, "r" ) ; 
if  (pif ile==NULL) 

{ 

printf ("\n  Error  opening  %s  for  reading.  Program  terminated.", 
inputf ilename) ; 

exit ( 1 ) ; 

} 

printf (" input  file  is  open  \n"); 
pofile=fopen (outputf ilename,  "a")  ; 
if (pof ile==NULL) 

{ 

printf ( "can ' t  open  %s  for  appending . \n" , outputf ilename)  ; 
exit (2 ) ; 

} 

printf ( "output  file  is  open  \n"); 

printf ("  Air  Volume=  %4.1f  mA3,  V=  %5.2f  m/s,  dTa/dday=  %5.2f  C/day 
\n" , Vps , V, Ta_rise  ); 

printf ("  wt=  %5.1f  kg,  ht=  %4.2f  m,  Bfat=  %4.1f  \n\n" , Wt, Ht, Bfat) ; 

{ 

fprintf (pofile, "\n\n  output  file  name:  %30s  \n",  outName) ; 
fprintf (pofile,  "  Air  Volume (mA3)=  %f,  V(m/s)=  %f,  dTa/dday=  %f 
\n" , Vps , V, Ta_rise  ); 

fprintf (pofile, "  wt(kg)=  %f,  ht (m) =  %f,  Bfat=  %f  \n\n" , Wt, Ht, Bfat) ; 
fprintf (pofile, "  hr  Ta  clo  RH  Tdp  Wo  Tc  Tsk  met 

sw  wet  HR  CWLos  CdripL  PSI  \n");  //header  for  ouput  file 
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fprintf (pofile, "  C  C  C  C 

g/min  bpm  L  L  \n"); 

} 

OutputHeader ( ) ;  //used  for  screen  display  only 

//TIM  Ta  clo  RH  Tdp  Wo  Tc  Tsk  met  sw  wet 

HR  CWLOS  PS I 

print f ( " \n%6 . If  %5.1f  %4.2f%4.0f  %5.1f  %5.3f  %5.2f  %5.2f  %4.1f  %5.2f 
%3.2f  %3 . Of  %8.2f%4.1f",TIM,Ta,clo,  RH, Tdpa, Wo, Tc, Tsk, met, Regsw/ 60  *Adu, wet, HR, 
CWLL,  PSIc)  ; 

fprintf (pofile, "\n%6. If  %5.1f  %4.2f%4.0f  %5.1f  %5.3f  %5.2f  %5.2f  %4.1f 
%5.2f  %3.2f  %3.0f  %8.2f  %8.2f  %4 . If " , TIM, Ta, clo, RH, Tdpa, Wo, Tc, Tsk, met, 
Regsw/60*Adu, wet, HR, CWLL, CdripL,  PSIc)  ; 

TIM  +=DT;  //  0+1  time  step  hr 

ReadTime+=DTM/ 60 ;  //time  for  next  read  in  hr  (default=  lhr) 

PrintTime+=DTMP/ 60 ;  //next  time  to  print  in  hours 
Tcl=Tsk;  //a  first  guess 

//MeasTime+=DTM* 60 ;  //  time  in  hr  for  next  measuredlnput 
time+=Duration;  //time  in  hr 

while  (TIM<=time) //  thermo-physiology  loop  start  measure  and  print  loop 

{ 

MTIM=TIM;  //TIM  in  hr  MTIM  in  hr 

if  (MTIM>=ReadTime) {  /*  if  current  time (MTIM) is  >  or  =  ReadTime 

it  reads  file.  if  not  it  assumea  conditions  are  unchanged*/ 
fscanf (pifile, "%101f%101f " , SExTime, SMetnew) ; 

ReadTime=ExTime+DTM/ 60 ; / /nextread  time  &  ExTime  in  hr 

} 

else  { 

Metnew=0 ; 

} 

if  (Metnew>0) 
met=Metnew; 

if  (Tsk>34  )  { 

clo=cloDisSub; 

} 

else  { 

clo=clos ; 

} 

FACL=1+ . 2*clo;  //estimate  of  outside  area  of  clothed  person 
Rcl= . 155*clo; 

Rpcl=  .153181*clo; 

/ /HR=HeartRate (met, Ta,  HRrest)  ; 

//met=MetT (HR, HRrest, Ta) ;  //from  Berglund  30th  ACEMB,1977 
Bf=  BreathRate (met, Wt, Adu) ;  //breathing  rate 

Bf c=Bf lowCore (Rm, HRRk) ;  //  ml/(min*mA2) 

Bfm=Bf lowMuscle (Rm, Bfc) ;  //  ml/ (min  mA2) 

SV=StrokeVolume (Adu, HRrest, Tsk) ;  //  ml 

Skbf=SkinBloodFlow (Tc, Tsk) ;  //  Liters/ (h  mA2) 

HR= (Skbf *1000/ 60+Bfc+Bfm) *Adu/ SV;  / /bpm 
HRRk=HR/HRkmax;  //fraction  of  maximum  heart  rate 
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/*  dry  and  evaporative  heat  transfer*/ 

Psk=SatVapPres (Tsk) ;  //vapor  press  of  water  on  skin  Torr 

he  =  Convection (met, V) ;  //  convection 

he=2.2*hc;  //  evaporation  watts/ (mA2  Torr) 

hr=4* . 725* (5 . 67E-08) *pow ( ( (Tcl+MRT) /2+273) , 3) ;  //  corrected  hr 
To=(hc*Ta  +  hr*MRT) / (hc+hr) ;  //operative  tempeature  C 

if  (code==0)  {  //used  intrinsic  clo 

//Fpcl=l/ ( 1+ . 15318 1 *he*FACL*clo) ;  //Berglund  1985,  for  IL=.45 
Rpbound=l/ (he*FACL) ; 

Rpclt=Rpcl+Rpbound; 

Rbound=l/ ( (hc+hr) *FACL) ; 

Rclt=Rbound+Rcl ; 

} 

else  {  //  used  total  clo 

Rclt=.155*Clot (code, V) ; 

Rpclt=VRes (code, V) ; 

} 

Emax= (Psk-Pa) /Rpclt;  / /watts/mA2 

Dry= (Tsk-To) / (Rclt) ;  //watts/mA2 

Tel  =  To  +Dry/ (FACE* (hr+hc) ) ;  //outer  clothing  temeperature  C 
/*  Thermal  Physiology  */ 

Res  =  Respiration (Rm, Ta, Pa, SWres) ;  //watts/mA2  Wres=g/ (h  mA2) 

Skbf =SkinBloodFlow (Tc, Tsk) ;  //Liters/ (h  mA2) 

alpha=Alpha (Skbf ) ;  //skin  fraction  of  total  weight 

Regsw=Sweat (Tc,  Tsk,  alpha)  ;  //  g/ (h  mA2) 

Ersw= . 68*Regsw;  //  watts/mA2  assuming  it  all  evaporates 

wet=Ersw/Emax; 

Edif  =  (1-wet) *. 06*Emax;  //watts/mA2  diffusion  through  dry  skin 
Wdif=Edif / . 68 ;  //  g/ (h  mA2)  rate  of  water  loss  by  diffusion 

if  (wet>=l)  { 
wet=l ; 

Wdif =0 . ; 

drip= (Wdif+Regsw) - (Emax/ . 68 ) ;  /*  g/ (h  mA2)  wasted  sweat 

driping  from  skin  */ 

Wsk  evap=Wdif +Regsw-drip;  /*  g/ (h  mA2)  water  evaporated  from 

skin* / 

} 

else  { 

drip=0 ; 

} 

//  skin  wettedness  from  sweating 

Wsk  evap=Regsw+Wdif ;  //water  vapor  loss  from  skin  g/ (hmA2) 

Esk= ( . 06+ . 94*wet) *Emax; 

RateOf TotalWaterLossPsqm=Regsw+Wres+Wdif ;  /*  g/ (h  mA2)  &  Wres  is 

g/ (h  mA2 )  */ 

RateOf TotalWaterLoss=RateOf TotalWaterLossPsqm*Adu;  /*  g/h 

includes  dripping*/ 
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Wevap=Regsw+Wdif +Wres ;  //g/ (h  mA2)  total  evaporation  rate  of 

water  from  person  per  mA2  */ 
PSIc=MoranPSI (  Tc,Tco,  HR,  HRrest); 

/ /PSIsk=BowlingPSI (Tsk, Tsko, HR,  HRrest)  ; 

HSCR  =  Core (Tc, Tsk, met, Res, SheatFlowCoreToSkin, Skbf, &Rm, &Bfat) ; 
HSSK  =  Skin (wet, Dry, Emax, SheatFlowCoreToSkin)  ; 

/*  thermal  capacity  */ 

TCCR=. 97* (1-alpha) *Wt; 

TCSK= . 97*alpha*Wt; 

if  (MTIM>=PrintTime) 

{ 

print f ( " \n%6  .  If  %5.1f  %4.2f%4.0f  %5.1f  %5.3f  %5.2f  %5.2f  %4.1f 
%5.2f  %3.2f  %3.0f  %8.2f  %4.1f",  " , TIM, Ta, clo, RH, Tdpa, Wo, Tc, Tsk, Rm/58 . 2 , 

Regsw/ 60*Adu, wet, HR, CWLL, PSIc) ; 

fprintf (pofile, "\n%6.1f  %5.1f  %4.2f%4.0f  %5.1f  %5.3f  %5.2f 
%5.2f  %4 . If  %5.2f  %3.2f  %3.0f  %8 . 2f %8 . 2f . %4 . If " , TIM, Ta, clo, RH, Tdpa, Wo, Tc, Tsk, 
Rm/ 58.2, Regsw/ 60*Adu, wet, HR, CWLL, CdripL, PSIc) ; 

PrintTime=ExTimetDTMP/ 60 ;  /*  next  print  time  &  ExTime  in  hr, 

DTMP  is  print  interval  in  hr  */ 

} 

TIM  +=DT;  //  hr 

//  Tc  and  Tsk  after  time  step  DT 
Tc  +=HSCR*Adu/TCCR*DT; 

Tsk  +=HSSK*Adu/TCSK*DT; 

Ta=Tat (Ta_rise/24) *DT; 

MRT=Ta; 

Cumu  1 at i veWa terLo s s Psqm=Cumul at i veWa terLo ssPsqm+RateOf Total Water Loss Psqm*DT; 

CumulativeWaterLoss=CumulativeWaterLoss+RateOf TotalWaterLoss*DT; 

//  g 

CWLL=CumulativeWaterLoss/1000;  //Liters  water  liq+evap 
CEvapWL  =  CEvapWL  +  Wevap*Adu*DT;  //g 

CEvapWLL  =  CEvapWL/1000;  //Liters  evap 

CdripL  =  CdripL  +  drip*Adu*DT/1000 ;  //Liters  drip 

II - 


Wo=Wi+CWLL/ma;  //  kgw/kga,  1  Lw  =  1  kgw 

Pa=Wo*Pbar* ( (Ta+273 ) / (Tao+273 ) ) / (Wot . 622 ) ;  /*  new  Pa  torr  with 

new  Pb=Pbar* (Ta+273) / (Tao+273) */ 


Psa=SatVapPres (Ta) ; 


if  (Pa>Psa) { 

Pa=Psa; 

Wos= . 622*Pa/ (Pbar-Pa) ;  //  kgw/kga  saturated  humidity  ratio 

Wo=Wos ; 


RH=Pa/Psa*100; 

Tdpa=4030 .183/ (18. 6686-log (Pa) ) -235; 
ExTime+=DT;  //hr 
}  //return  to  while 


fclose (pif ile) ;  //close  inputfile 
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fclose (pof ile) ;  //close  outputfile 
return  0  ; 

} 


void 

{ 

OutputHeader ( ) 

printf ( "\n  hr  Ta 

clo 

RH  Tdp  Wo  Tc 

Tsk 

met  sw 

wet 

HR 

bpm 

CWLoss  PSI  \n" )  ; 
printf ( "  C 

L\n" ) ; 

C  kgw/kga  C 

C 

g/min 

return  ; 

} 


double  SatVapPres (double  T) 

{ 

double  Pst; 

Pst  =  exp (18 . 6686- (4030 . 183/ (T+235 .))) ;  //mmhg 
return  Pst; 

} 

double  Convection (double  met, double  V) 

{ 

double  CHCA,  CHCV,  CHCmin,  he; 

CHCmin  =  3.0; 
if  (met>=l . 1 )  { 

CHCA  =  5 . 66*pow ( (met  -  0.85),  0.39);  /*  he  due  to  activity  */ 

} 

else  { 

CHCA  =  5.66*pow((l.  -  0.85),  0.39);  /*  he  due  to  activity  */ 

} 

CHCV  =  8 . 6*pow (V, 0 . 53 ) ;  /*  he  due  to  air  speed  V  in  m/s  */ 

if  (CHCV  >=  CHCmin) 


} 


else 

(CHCV  =  CHCmin) ; 
if  (CHCV  >=  CHCA) 
he  =  CHCV; 

else 


he  =  CHCA; 

return  he; 


double  Respiration (double  Rm,  double  Ta,  double  Pa, double*Wres) 

{ 

double  Res , Eres , Cres ; 

Eres  =  0 . 0023*Rm* (44 . -Pa) ;  /*  watts/m2  */ 

Cres  =  0 . 0014*Rm* (34 . -Ta) ;  /*  watts/m2  */ 

*Wres=Eres/ . 68 ;  /*  g/ (h  mA2)  */ 

Res  =  Eres  +  Cres; 
return  Res; 

} 


double  DuBoisSkinArea (double  Wt,  double  Ht) 

{ 

double  Adu;  /*  mA2  */ 

Adu=0 . 202*pow (Wt, 0 . 425) *pow (Ht, 0 . 725) ; /*Wt  weight  in  kg,  Ht  height  in 

m.  */ 
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return  Adu; 

} 

double  BflowCore (double  Rm,  double  HRRk) 

{ 

double  BFcr , BFcrmax; 

BFcrmax= . 92*Clave;  //  ml/(mA2*min) 
if (Rm/58 . 2<=1 . 3) {  //if  at  rest 
if  (HRRk> . 25 ) 

BFcr=BFcrmax* (1-1 .266* (HRRk- . 25 ) ) ; 

else 

BFcr=BFcrmax; 

} 

else  //during  exercise 

{ 

if  (HRRk>. 39) 

BFcr=BFcrmax* (1-1 .0866* (HRRk- . 39)  )  ; 

else 

BFcr=BFcrmax; 

} 

return  BFcr;  //ml/ (min*mA2) 

} 

double  BflowMuscle (double  Rm, double  Bfc) 

{ 

double  mcr, mmu, BFmu; 

mcr=Bf c*BL02cap*F02ExtCr/ 02equiv; 

mmu=Rm-mcr ; 

BFmu=mmu*02equiv/ (F02ExtMu*BL02cap) ;  /*muscle  bloodflow  in 

mlblood/ (min*mA2)  */ 

return  BFmu; 

} 

double  StrokeVolume (double  Adu,  double  HRrest, double  Tsk) 

{ 

double  SV1 , SV; 

SVl=Clave*Adu/HRrest;  //stroke  volume  ml 
if  (Tsk<=33) 

SV=SV1 ; 
else 
{ 

if  (Tsk<38 ) 

SV=SVl-5* (SV1-85) * (Tsk-33) /45; 
else 

SV=SVl-5* (SV1-85) * (38-33) /45; 

} 

return  SV; 

} 

double  Shiver (double  Tc,  double  Tsk, double*Bfat) /*Tikusis  &  Stolwijk  Models*/ 

{ 

double  shiver=0;  //  Bfat  is  %  body  fat 

if  (Tc<TTCR  | |  Tsk<33)  /*  does  shiving  occur?  Tc<37/*/ 

{  if  (Tc<TTCR  &&  Tsk<33)  /*shiver  diven  by  core  and  skin*/ 

shiver= (156* (TTCR-Tc) +47* (33-Tsk) -1 . 57*pow ((33- 
Tsk) , 2) ) /pow (*Bfat,  0 . 5) ;  /*Tikusis  and  Giesbrecht,  1999,  15%BFat  */ 
else 
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if (Tc<TTCR)  /*shiver  driven  core  only*/ 
shiver= (156* (TTCR-Tc) ) /pow (*Bfat, 0 . 5) ; 
else 

shiver= (47* (33-Tsk) -1. 57*pow ( (33-Tsk) ,2) ) /pow (*Bfat, 0 . 5) ; 

//driven  by  skin  only 

} 

//shiver=0;  /*Use  to  disable  shiver*/ 

return  shiver; 

} 

double  MetT (double  HR,  double  HRrest, double  T)  / /Berglund, 30th  ACEMB  1977 

{ 

double  met,HRratio;  /*  applicable : 20<=T<=40,  0.6clo,  v~=1.25m/s, 

1 ,2<+HRratio<=2 . 1, Tdp<=20C*/ 

HRratio=HR/HRrest; 
if  (HRratio>=l . 2 )  { 

if  (T>=20 ) 

met=0 . 68+4 . 69* (HRratio  -1)-  0 . 052* (HRratio-1 ) * (T-20 ) ; 
else 

met=0 . 68+4 . 69* (HRratio  -1); 

} 

else  { 

met=HRratio; 

} 

return  met; 

} 

double  BreathRate (double  Xmet,  double  Wt,  double  Adu)  /*aerobic  breathing 
rate  for  metabolism*/ 

{ 

double  Bf,  V02pkg; 

V02pkg=Xmet*58 . 2*02equiv*Adu/Wt; 

Bf= . 3957*V02pkg+18 .735; 
return  Bf; 

} 

double  HeartRate (double  met, double  T,  double  HRrest)  /*from  Berglund  30th 

ACEMB, 1977  */ 

{ 

double  HRratio, HR; 

HRratio=l+ (met- . 68 ) / (4.69-  0.052*  (T-20) ) ; 

HR=HRrest*HRratio; 
return  HR; 

} 

double  Alpha (double  Skbf ) 

{ 

double  alpha; 

alpha=0.04177+. 74518/ (Skbf+0 . 585417) ; 
return  alpha; 

} 

double  Sweat (double  Tc,  double  Tsk,  double  alpha) 

{ 

double  regsw,  Tmb; 
regsw=0 ; 

Tmb= (1-alpha) *Tc  +  alpha*Tsk; 
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if  ( (Tmb>TTBM) && (Tsk>TTSK) ) 

regsw=CSW* (Tmb-TTBM) *exp ( (Tsk-TTSK) /10 . 7) ; 
else  if  ( (Tmb>TTBM) && (Tsk<=TTSK) ) 
regsw=CSW* (Tmb-TTBM) ; 

if  (regsw>667)  //  max  sweat  rate  limit 

regsw=667;  //  regsw  max=667g/ (h  mA2)=ll.lg/ (min  mA2) !=20g/ (min  mA2) 
return  regsw;  //g/(h  mA2) 

} 

double  SkinBloodFlow (double  Tc,  double  Tsk) 

{ 

double  Colds=0; 
double  Skbf,  WarmC=0; 
if  (TskCTTSK) 

Colds=TTSK-Tsk; 
if  (Tc>TTCR) 

WarmC=Tc-TTCR; 

Skbf= (SKBFN+CDIL*WarmC) / (l+CSTR*Colds)  ;  //  Liters/ (h  mA2 ) 

if  ( Skbf >Skbfmax) 

Skbf=  Skbfmax; 
if  (Skbf<Skbfmin) 

Skbf=  Skbfmin; 

return  Skbf;  //L/ (h  mA2) 

} 

double  Core (double  Tc, double  Tsk, double  met, double  Res, 

double*heatFlowCoreToSkin,  double  Skbf , double*Rm, double*Bfat) 

{ 

double  Rmet,  Hfcrsk, HSCR, shiver ; 

shiver=Shiver (Tc,  Tsk,  Bf at)  ;  /*  watts/mA2  */ 

Hfcrsk= (CMIN+CB*Skbf ) * (Tc-Tsk) ;  /*  watts/mA2  */ 

Rmet  =  58.2*met  +  shiver;  /*  metabolic  heat  produced  watts/mA2  */ 
HSCR=Rmet-Hf crsk-Res ;  /*  rate  of  heat  storage  in  core  watts/mA2  */ 

*heatFlowCoreToSkin=Hf crsk; 

*Rm=Rmet; 

return  HSCR;  //rate  of  heat  storage  in  core  compartment  w/mA2 

} 

double  Skin (double  wet, double  Dry, double  Emax, double*heatFlowCoreToSkin  ) 

{ 

double  Esw, Edif f , Esk, HSSK; 

Esw=wet*Emax; 

Edif f = . 0  6* ( 1-wet ) *Emax; 

Esk=Esw+Edif f  ; 

HSSK=*heatFlowCoreToSkin-Dry-Esk; / *rate  of  heat  storage  in  skin  w/mA2*/ 
return  HSSK; 

} 

double  MoranPSI (double  Tc, double  Tco, double  HR,  double  HRrest) 

{ 

double  PSI; 

PS 1=5* ( (Tc-Tco) / (39 . 5-Tco) + (HR-HRrest) / (180-HRrest) ) ; 
return  PSI; 

} 

/*  This  is  the  End  */ 
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/* 

Example  of  Input  and  resulting  output  files  of  Sub  ICDA 


Input  file  name:  Subin3.txt 


hr 

met 

l . 

l . 

2  . 

l . 

3. 

l . 

4  . 

l . 

5. 

l . 

6. 

l . 

7  . 

l . 

8  . 

l . 

Output : 

output  file  name:  Sub0ut3bVps  RH.txt 

Air  Volume(mA3)=  360.000000,  V(m/s)=  0.050000,  dTa/dday=  1.110000 
wt(kg)=  82.000000,  ht (m) =  1.770000,  Bfat=  15.000000 


hr 

Ta 

C 

clo 

RH 

Tdp 

C 

Wo 

Tc 

C 

Tsk 

C 

met 

sw  wet 
g/min 

HR 

bpm 

CWLos 

L 

CdripL 

L 

PSI 

0.0 

24.0 

0.20 

60 

15.8 

0.011 

36.90 

33.00 

1.0 

0.00 

0.06 

70 

0.0 

0.0 

0.0 

6.0 

24.3 

0.20 

62 

16.5 

0.012 

36.90 

32.48 

1 . 1 

0.00 

0.00 

74 

0.2 

0.0 

0.2 

12.0 

24 . 6 

0.20 

63 

17.1 

0.012 

36.89 

32.59 

1.1 

0.00 

0.00 

74 

0.4 

0.0 

0.2 

18.0 

24 . 8 

0.20 

65 

17 . 8 

0.013 

36.89 

32.70 

1.1 

0.00 

0.00 

74 

0 . 6 

0.0 

0 . 1 

24.0 

25.1 

0.20 

66 

18.3 

0.013 

36.88 

32 .81 

1.0 

0.00 

0.00 

73 

0 . 8 

0.0 

0 . 1 

30.0 

25.4 

0.20 

67 

18 . 9 

0.014 

36.88 

32 . 92 

1.0 

0.00 

0.00 

73 

1.0 

0.0 

0 . 1 

36.0 

25.7 

0.20 

68 

19.4 

0.014 

36.87 

33.06 

1.0 

0.00 

0.00 

73 

1.2 

0.0 

0 . 1 

42.0 

25.9 

0.20 

70 

19.9 

0.015 

36.88 

33.34 

1.0 

0.00 

0.00 

73 

1.4 

0.0 

0 . 1 

48.0 

26.2 

0.20 

71 

20.4 

0.015 

36.88 

33.57 

1.0 

0.03 

0.00 

74 

1.5 

0.0 

0 . 1 

54.0 

26.5 

0.20 

72 

21.0 

0.015 

36.88 

33.63 

1.0 

0.11 

0.02 

74 

1 . 8 

0.0 

0 . 1 

60.0 

26.8 

0.20 

73 

21.5 

0.016 

36.88 

33.70 

1.0 

0 .19 

0.03 

74 

2.0 

0.0 

0 . 1 

66.0 

27 . 1 

0.20 

75 

22 . 1 

0.017 

36.88 

33.75 

1.0 

0.28 

0.05 

74 

2.2 

0.0 

0 . 1 

72.0 

27.3 

0.20 

76 

22 . 8 

0.017 

36.89 

33.81 

1.0 

0.37 

0.07 

74 

2.5 

0.0 

0.2 

78.0 

27 . 6 

0.20 

78 

23.5 

0.018 

36.89 

33.86 

1 . 0 

0.46 

0.09 

74 

2 . 8 

0.0 

0.2 

84.0 

27 . 9 

0.20 

80 

24.2 

0.019 

36.90 

33.91 

1.0 

0.55 

0 . 12 

74 

3.1 

0.0 

0.2 

90.0 

28.2 

0.20 

82 

24 . 9 

0.020 

36.91 

33.97 

1.0 

0.65 

0.14 

74 

3.5 

0.0 

0.2 

96.0 

28.4 

0.20 

85 

25.6 

0.021 

36.91 

34.03 

1.0 

0.75 

0 . 18 

74 

3.9 

0.0 

0.2 

102.0 

28.7 

0.20 

87 

26.4 

0.021 

36.92 

34.08 

1.0 

0 . 84 

0.21 

74 

4.3 

0.0 

0.2 

108.0 

29.0 

0.20 

90 

27.2 

0.023 

36.93 

34 . 14 

1.0 

0 . 94 

0.26 

75 

4.7 

0.0 

0.3 

114.0 

29.3 

0.20 

93 

27 . 9 

0.024 

36.93 

34 .19 

1.0 

1.05 

0.31 

75 

5.2 

0.0 

0.3 

120.0 

29.5 

0.20 

95 

28.7 

0.025 

36.94 

34.25 

1.0 

1 . 15 

0.38 

75 

5.6 

0.0 

0.3 

126.0 

29.8 

0.20 

98 

29.5 

0.026 

36.95 

34.31 

1.0 

1.25 

0.47 

75 

6.1 

0.0 

0.3 

132.0 

30.1 

0.20 

100 

30.1 

0.027 

36.96 

34.36 

1.0 

1.35 

0.56 

75 

6.7 

0.0 

0.3 

138.0 

30.4 

0.20 

100 

30.4 

0.028 

36.96 

34.41 

1.0 

1.44 

0 . 63 

75 

7.2 

0.0 

0.4 

143.9 

30.7 

0.20 

100 

30.7 

0.028 

36.97 

34.46 

1.0 

1.52 

0.70 

75 

7 . 8 

0.0 

0.4 

149.9 

30 . 9 

0.20 

100 

30 . 9 

0.029 

36.98 

34.50 

1.0 

1 . 61 

0.78 

75 

8.4 

0.0 

0.4 

155.9 

31.2 

0.20 

100 

31.2 

0.029 

36.98 

34.55 

1.0 

1.70 

0 . 87 

76 

9.0 

0.0 

0.4 

161 . 9 

31.5 

0.20 

100 

31.5 

0.030 

36.99 

34 . 60 

1.0 

1.79 

0 . 98 

76 

9.7 

0.0 

0.4 

167 . 9 

31 . 8 

0.20 

100 

31 . 8 

0.030 

37.02 

34 . 80 

1.0 

2 .16 

1.00 

76 

10.4 

0.0 

0.5 

173.9 

32.0 

0.20 

100 

32.0 

0.031 

37.07 

35.05 

1.0 

2 . 67 

1.00 

77 

11.3 

0.2 

0 . 6 

179.8 

32.3 

0.20 

100 

32.3 

0.031 

37 . 12 

35.30 

1.0 

3.24 

1.00 

78 

12.3 

0 . 6 

0 . 8 

185.8 

32 . 6 

0.20 

100 

32 . 6 

0.032 

37 .18 

35.56 

1.0 

3.87 

1.00 

79 

13.6 

1.2 

0 . 9 

191 . 8 

32 . 9 

0.20 

100 

32 . 9 

0.032 

37.26 

35.81 

1.0 

4.59 

1.00 

80 

15.2 

2 . 1 

1 . 1 

197 . 8 

33.1 

0.20 

100 

33.1 

0.033 

37.34 

36.06 

1.0 

5.42 

1.00 

81 

17.0 

3.2 

1.4 

203.8 

33.4 

0.20 

100 

33.4 

0.033 

37.44 

36.31 

1.0 

6.37 

1.00 

83 

19.1 

4.7 

1 . 6 

209.8 

33.7 

0.20 

100 

33.7 

0.034 

37.56 

36.56 

1.0 

7.45 

1.00 

84 

21 . 6 

6.5 

1 . 9 
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215. 

.7 

34  . 

.0 

0  . 

.20 

100 

34  . 

.0 

0  . 

.034 

37  . 

.  69 

221 . 

.7 

34  . 

.3 

0  . 

.20 

100 

34  . 

.3 

0  . 

.035 

37  . 

.  84 

227  . 

.7 

34  . 

.5 

0  . 

.20 

100 

34  . 

.5 

0  . 

.036 

38  . 

.00 

233. 

.7 

34  . 

.  8 

0  , 

.20 

100 

34  . 

.  8 

0  . 

.036 

38  . 

.  18 

239. 

.7 

35. 

.1 

0  . 

.20 

100 

35. 

.1 

0  . 

.037 

38  . 

.37 

245. 

.7 

35. 

.4 

0  . 

.20 

100 

35. 

.4 

0  . 

.037 

38  . 

.58 

251 . 

.  6 

35. 

.  6 

0  . 

.20 

100 

35. 

.  6 

0  . 

.038 

38  . 

.  84 

257  . 

.  6 

35. 

.  9 

0  . 

.20 

100 

35. 

.  9 

0  . 

.039 

39. 

.  10 

263. 

.  6 

36. 

.2 

0  . 

.20 

100 

36. 

.2 

0  . 

.039 

39. 

.35 

269. 

.  6 

36. 

.5 

0  . 

.20 

100 

36. 

.5 

0  . 

.040 

39. 

.  61 

275. 

.  6 

36. 

.7 

0  , 

.20 

100 

36. 

.7 

0  . 

.040 

39. 

.  87 

281 . 

.  6 

37  . 

.0 

0  . 

.20 

100 

37  . 

.0 

0  . 

.041 

40  . 

.  13 

287  . 

.  6 

37  . 

.3 

0  . 

.20 

100 

37  . 

.3 

0  . 

.042 

40  . 

.39 

293. 

.  6 

37  . 

.  6 

0  . 

.20 

100 

37  . 

.  6 

0  . 

.042 

40  . 

.  65 

299. 

.  6 

37  . 

.  9 

0  . 

.20 

100 

37  . 

.  9 

0  . 

.043 

40  . 

.  92 

305. 

.  6 

38  . 

.1 

0  . 

.20 

100 

38  . 

.  1 

0  . 

.044 

41 . 

.  18 

311 . 

.  6 

38  . 

.4 

0  . 

.20 

100 

38  . 

.4 

0  . 

.045 

41 . 

.44 

317  . 

.  6 
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